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SLEW MANEUVERS OF SPACECRAFT CONTROL 
LABORATORY EXPERIMENT (SCOLE) 

Y. P. Kakad 

Dept, of Electrical Engineering 
University of North Carolina at Charlotte 
Charlotte, NC 28223 


SUMMARY 

This is the final report on the dynamics and control of slew maneuvers of the 
Spacecraft COntrol Laboratory Experiment (SCOLE) test facility initiated and 
developed by the Spacecraft Controls Branch at NASA Langley Research Center. 
The report includes the work done under grant period from November 1st., 1984 
througn April 30th, 1992. Although the final experimental work to be performed 
on the test facility is scheduled to be done in May 1992, the report includes all 
theoretical and numerical work upto the closing date of the grant. A supplemen- 
tary report will be issued on the experimental work after the remaining work is 
completed. In the appendix, two current publications on the grant work are 
included and the other publications based on the grant work are cited in the refer- 
ences. 

The report documents the basic dynamical equation derivations for an arbi- 
trary large angle slew maneuver as well as the basic decentralized slew maneuver 
control algorithm. 

The dynamics and control of slew maneuver of NASA Spacecraft COntrol 
Laboratory Experiment (SCOLE) test facility are developed in terms of an arbi- 
trary maneuver about any given axis. The set of dynamical equations incorporate 
rigid-body slew maneuver and three-dimensional vibrations of the complete 
assembly comprising the rigid shuttle, the flexible beam, and the reflector with an 
offset mass. The analysis also includes kinematic nonlinearities of the entire 
assembly during the maneuver and the dynamics of the interaction between the 
rigid shuttle and the flexible appendage. The final set of dynamical equations 
obtained for slewing maneuvers are highly nonlinear and coupled in terms of the 
flexible modes and the rigid-body modes. 

The equations are further simplified and evaluated numerically to include the 
first ten flexible modes and the SCOLE data to yield a model for designing con- 
trol systems to perform slew maneuvers. 

The control problem formulation incorporates the nonlinear dynamical equa- 
tions and is expressed in terms of a two-pomt boundary value problem utilizing a 
quadratic type of performance index. 

The two-point boundary value problem is solved as hierarchical control 
problem with the overall system being split in terms of two subsystems, namely 
the slewing of the entire assembly and the vibration suppression of the flexible 
antenna. The coupling variables between the two dynamical subsystems are 
identified and these two subsystems for control purposes are treated 
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independently in parallel at the first level. Then the state- space trajectory of the 
combined problem is optimized at the second level. 


1. INTRODUCTION 

The primary control objective of the Spacecraft Control Laboratory Experi- 
ment (SCOLE) is to direct the RE Line-Of- Sight (LOS) of the antenna-like 
configuration towards a fixed target under the conditions of minimum time and 
limited control authority [1]. This problem of directing the LOS of antenna- like 
configuration involves both the slewing maneuver of the entire assembly and the 
vibration suppression of the flexible antenna-like beam. The study of ordinary 
rigid-body slew maneuvers has received considerable attention in the literature 
[2,3] due to the fact that any arbitrary large-angle slew maneuver involves 
kinematic nonlinearities. This is further complicated in the case of SCOLE by vir- 
tue of a flexible appendage deployed from the rigid space shuttle. The dynamics 
of arbitrary large-angle slew maneuvers of SCOLE model are derived in this 
report as a set of coupled equations with the rigid-body motions including the 
nonlinear kinematics and the vibratory equations of the flexible appendage. 

The dynamical equations of slewing maneuvers of this large flexible space- 
craft are developed by writing the total kinetic and potential energy expressions 
for the entire system. The energy expressions are further utilized in formulating 
Lagrange’s equations which are expressed in terms of non-generalized co- 
ordinates using an inertial co-ordinate system and a body-fixed co-ordinate sys- 
tem at the point of attachment of the flexible beam to the shuttle. The generic 
model used for this analysis consists of a distributed parameter beam with two 
end masses. The three dimensional linear vibration analysis of this free-free beam 
model with end masses [4] is incorporated together with rigid-slewing maneuver 
dynamics which are written in terms of four Euler parameters [5] and angular 
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rotation about an arbitrary axis of rotation to yield the final set of highly nonlinear 
and coupled equations. In the derivation of the equations, it is assumed that the 
vibratory analysis is for small motions. These nonlinear and coupled dynamical 
equations are used in this article to study the slew maneuver control in terms of a 
hierarchical feedback control scheme. 

The control problem of slewing maneuvers of this large flexible spacecraft is 
developed by using the two-point boundary value problem in terms of the rigid- 
body slewing and the vibration suppression of the flexible appendage as two 
separate dynamical subsystems. A decentralized optimal control scheme is util- 
ized in order to solve individual boundary-value problem for each of the two sub- 
systems by defining their state variable models and incorporating the coupling 
variables between the two subsystems in these models. Also, the boundary condi- 
tions of the overall system are reworked in terms of boundary conditions of each 
subsystem. A quadratic performance index is utilized for the overall system and is 
subsequently expressed in terms of a sum of two individual performance indices 
of the subsystems. 

The basic algorithm for obtaining an optimal closed-loop state feedback 
sceme involves using a trajectory in terms of a vector of Lagrange multipliers as 
an initial guess at level two. This is used at level one in quasilinearization appli- 
cation. 

The two-point boundary value problem for each subsystem is solved at level 
one by using a quasilinearization technique as a trajectory optimization problem. 
In the quasilinearization procedure, starting from an initial guessed state trajec- 
tory, successive linearizations are performed of state equations in such a way that 
the final solution of the state trajectory is within an acceptable degree subject to 
boundary conditions. The state vector definition at this level is an augmented 
state vector which includes both system states and costates. 
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These optimum solutions of the system trajectories are utilized at level two 
to yield the updated trajectory of the vector of Lagrange multipliers of the overall 
system to be used for quasilinearization process at level one. The basic steps of 
the algorithm are repeated to optimize this second level trajectory with respect to 
prespecified error criterion to obtain an optimal feedback law. 

2 . LIST OF SYMBOLS 


a (z) Position vector of mass element on the beam from the point 

of attachment 

B Damping matrix 

C Inertial frame to body-fixed frame transformation 

c Position vector from the point of attachment to the mass center 

of the beam 

D Mass density of the beam 

d(z,t) Displacement vector of mass element in the body-fixed frame 
E Modulus of Elasticity 

Fo(t) Force applied at the orbiter mass center 
F 2(0 Force applied at the reflector mass center 
G 0 (t) Moment applied about the orbiter mass center 
Modulus of rigidity for the beam 
I Beam cross section moment of inertia 

l x Beam cross section moment of inertia, roll bending 

I y Beam cross section moment of inertia, pitch bending 

I 0 Equivalent mass moment of inertia 

/ 1 Mass moment of inertia matrix of the shuttle 

1 2 Mass moment of inertia matrix of the reflector 

J Mass moment of inertia matrix of the beam 

K Stiffness matrix 

L The Length of the beam 
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M Angular velocity vector transformation 

My Effective moment applied at the reflector c.g. 

m Total mass of the flexible beam 
m i Mass of the orbiter 

m2 Mass of the reflector 

N The total number of subsystems 

n The maximum number of modes considered 
Q The generalized force vector 

qi Generalized coordinates 

R Position vector of the mass center of the orbiter in the inertial 
~ frame 

r Position vector from the orbiter mass center to the point of 

attachment 

r x x co-ordinate of the reflector mass center in the body-fixed 

frame 

r v y co-ordinate of the reflector mass center in the body-fixed 

frame 

T Total Kinetic Energy 

U Total Potential Energy 

u x (z,t ) The beam deflection in x direction referred to the body- 

fixed frame 

u y (z,t ) The beam deflection in y direction referred to the body- 

fixed frame 

My(z,r) The torsional deflection about z axis in the body- 

fixed frame 

Uj Control vector of ith system 

V Velocity vector of the mass center of the orbiter in the body- 

fixed frame 

V 0 Velocity vector of the point of attachment in the body- 

fixed frame 

Xj State vector of ith system 

Z( Vector of interconnecting variables 
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p Mass per unit length of the flexible beam 

Y Vector representing the axis rotation during the slew 

- maneuver 

i th Eigenfunction corresponding to u x 
<|y i th Eigenfunction corresponding to u y 
<|> v ,- i th Eigenfunction corrsponding to 
0 The attitude of the orbiter in the inertial frame 
£ Slew Angle 

0 ) The angular velocity of the orbiter in the inertial 

frame 

Q The angular velocity of the reflector in the inertial 

— frame 

e Vector of Euler parameters 

8(z— z j) Direc delta function 

<b(A.) Dual functional for two-point boundary value problem 

X Vector of Lagrange multipliers 

£ Damping ratio 
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3 . ANALYTICS 

Co-ordinate Systems 

The motion of SCOLE assembly when considered as a rigid body in space 
has six dynamic degrees of freedom: three of these define the location of the mass 
center, and three define the orientation (attitude) of the body. The motion of this 
rigid body is goverened by newtonion laws of motion expressed in terms of 
changes in linear momentum and angular momentum. These relationships are 
valid only when the axes along which the motion is resolved are an inertial frame 
of reference [ 9 , 10 ]. To define the orientation of the orbiter in space, a set of 
orthogonal axes fixed in the body is utilized. Then the attitude of the orbiter is 
defined in terms of the angles (0 1,62163) between the body- fixed axes and the 
inertial co-ordinate axes. The body-fixed frame origin is located at the point of 
attachment of the flexible appendage with the rigid shuttle for this analysis (Fig. 
1 ). 

The transformation from the inertial frame to the body-fixed frame is given 
by the matrix, C as developed in figure 2 where if "it j» represent the dexteral 
set of orthogonal unit vectors fixed in the body- fixed frame and 0i is the rotation 
about i, 62 is the rotation about j and 63 is the rotation about k. These rotations 
are carried out successively as shown in figure 1 and the matrix C is given as 



COS03 

sin 03 0 


cos 02 0 -sin 02 


1 0 0 

c = 

-sin 03 cos 03 0 


0 1 0 


0 cos0i sin0i 


0 

0 1 


sin 02 0 cos 02 


0 -sin 0i cos0i 


Thus C T is obtained as 


C r = 


cos02cos03 -cos02sin03 sin02 

sin0isin02cos03+sin03cos0i -sin0isin02sin03+cos03cos0i -sin0icos02 
-cos0i sin02cos03 +sin03sin0i cos0i sin02sin03 +cos03sin0i cos0icos02 


( 2 ) 
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In order to completely define the attitude (orientation), it is needed to relate 
the rotation angles 0i, 02 , and 03 to the angular velocity components (coi, 0 ) 2 , 0)3) 
of the orbiter. One way of obtaining the required relations is via body-three 
angles method [5] which was utilized in developing C matrix in equation (1) and 
these relations are 


01 = ( COiCOS 03 - (O2sin03)/cos02 

02 = ( (OiSin 03 +CO 2 COS 03 ) (3) 

03 = (-coicos 03 +co 2 sin 03 )tan 02 +o )3 

Thus, the angular velocity of the orbiter can be obtained in the inertial frame 
by means of the following transformation 


to = M r 0 

where the transformation M r is given as 


(4) 


M t = 


COS 02 COS 03 sin 03 0 
-cos 02 sin 03 cos 03 0 
sin 02 0 1 


(5) 


Although the body-three angles method is used here for obtaining the 


transformations C and M, there are three other methods which can be used to 


obtain the same transformations. A detailed discussion of all the methods is given 
in reference [5] and a summary of the transformations using the remaining three 
methods is given in the Appendix. 


Kinetic Energy 


If the position vector of the mass center of the orbiter in the inertial frame 
(Fig. 3), R , is given as 
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R = 


Rx 

Ry 

Rz 


then the velocity of the mass center in the inertial frame is 


( 6 ) 



Rx 

Ry 

Rz 


This velocity can be transformed in the body-fixed frame as 


(7) 


m = c 


Rx 

Ry 

Rz 


( 8 ) 


The velocity of the point of attachment in the body-fixed frame is 


V 0 = K+coxr (9) 

where r is the vector from orbiter mass center to the point of attachment. 

Defining the position vector (Fig. 4), a, of a mass element on the beam from 
the point of attachment (origin of the body-fixed frame) before deformation as 


a = 


F01 

0 

z 


( 10 ) 


and the displacement vector of this mass element as 


d(z,t) = 


u x (z,t) 

Uy(Z,t) 

0 


(ID 


the position vector after deflection is given as a+d. The kinetic energy in the 



- 10 - 


beam [6] is 

T i = (1/2) mVlVo + (1/2) (0 T [J](0-mVo T [c](a+(V2) J d T ddm 

+ yl J ddm+u/ J addm+( 1/2) J ^ ii x ' ii y ' dl 

where the vector c is from the point of attachment to the mass center of the beam 
and if it is assumed that the beam is a thin rod, then it is given as 


Mjc 

Uy' 

Wy 


( 12 ) 


c = 



= (1/m) J adm = 


0 ' 
0 

-LI 2 


(13) 


and using the skew symmetric form for the vector cross product for any two vec- 
tors c and to (in the same reference frame) as 

cx co = [ c ] co 


c = 


0 -c 


Cy C 


z y 

o -Cx I 
0 


also, the moment of inertia matrix is given as 


7 = (l/3)pL 3 


10 0 ' 
0 1 0 
0 0 0 


(14) 


(15) 


where p is the mass per unit length of the beam. The last term in the equation (12) 
corresponding to torsional motion is given as 
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l/2(p6ts , )5- 2 0 0 


u x ' 

0 l/2(pds)s 2 0 


Uy ' 

0 0 0 


J 

■ 


M\j/ 


( 16 ) 


The kinetic energy equation (12) can be simplified as 


T i 


( 1 / 2 ) pLV%V 0 +(l/6) pL 3 f to? -Koi 1 -pLVZca+pL £ qj 

1 J i=i 


+Vjd+co r p+(l/ 4 ) p 


ft 2 n 2 

£P5i<?i + 2>6i<?j 
i=l i =1 


( 17 ) 


where 


n 

u x = ^ §xi( s )QiO ) 
i =1 
n 

Uy = y ^>vi(<y)<7i(/) 
i=l 

u x ~ tyxi ( s )tfiO ) 

1=1 

My = (^)^i(O 

1=1 
n 

<=1 

L 

P \i = \$xi(s)ds 
0 
L 

P2i = l$yi(s)ds 
0 
L 

Pli = 

0 

L 

P4i = js$ y i(s)ds ( 18 ) 


0 
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and 


a(0 


P 


■ n 

i = 1 
n 

'LPHQi 

1=1 

0 


( 19 ) 


m = p 


■ n 

'LPMRi 

1=1 

n 

EP3i<7j 
i =1 

0 


( 20 ) 


The expressions for p u, pn> Pn, P 4 i, P 5 i , and psi are developed as follows. Note 
that 


(j >xi(s) = Axjsinfiis+Bxicosfiis+Cxismhfiis+Dxicoshfiis 


where P, = 



1/4 


El 


Since for SCOLE configuration EI X = EI y and P iJC = P,y , El and P, are used 
for both <t>ja(j) and (Jyfa ). However, this may not be true for other configurations. 


Pli = 


1 



L 

p i i = J^(*y)^ 

0 

L 

Pli = J<|> 7 /(s)ds 
0 

| -AjaCOsPiL +B X i sinaj +Ca cosh P t L 

+Dxi sinhP,L + Axi —C# 


(21 A) 


Defining otj = P;L 
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Pn = ~~ [ -A^cosa i +5^sina J +C JC iCosha,+Z) A3 sinha,+A^-C XJ ] 
Ctj/> L J 

similarly. 


P2i = 


PH = 


2 J" A^cosp/L +5y,sinP J L 
P/L L 

+C-y t coshPjL +Dyi sinhP/L +A y ~ C y i j 
£ -A^cosa,- +B > ,sina J +Cy/Cosha, +D yi sinha, +A y /-Cy/ j 


Pli - 


o 

L 


Pm = is$yi(s)ds 
0 


and these can be given as 

sinp/L LcosP,L 


/?3/ = A^ 


P3i=A 


xi 


xi 


P f Pi ' +B 

Lcoshp/L sinhp,L 

;* a- ~ p? 

L 2 sina,- L 2 cosa, 


of 


a i 


' XI 


cosP,L LsinP/L i 


P t 


+D 


XI 


+B 


xi 


ft p? 

Lsinhp,L coshP/L _1_ 

p. ' p? + p? 


L 2 cosa, L 2 sina, £ 2 
•+ 


a 1 


J _ 

L 2 cosha,- L 2 sinha/ 

k 


+ Dxi 

«i af 



«i a 2 

L 2 sinha/ L 2 cosha, L 


a i 


a ? 


a 


Similarly, 


PM - Ayi 


sinP,L LcosP/L , 
Lcoshp,L sinp/L 


'j* 


PAi — Ayi 


Pi 


P? 


L 2 sina,- L 2 cosa/ 


a} 


a, 


cosp/L L sinp/L i 
+- 


Pf 


+L> 




+5 




p/ p? 

LsinhP/L coshp/L l 

ft " P, ? + Pi 


L 2 cosa, L 2 sina,- /, 2 
a 2 a/ a 2 


(21B) 


(22A) 

(22B) 


(23A) 


(23B) 


(24A) 
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'yi 


L 2 cosha z L 2 sinha/ 


a. 


a 1 


+D 


yi 


L 2 sinha z L 2 cosha, jJ 1 


P Si 


L 

= 1 [sQxiV) j 

L 

P6i = J S^yiXs) 
0 1 


a i 
2 

ds 

2 

ds 


2 + 2 
a f a t 


and these can be shown to be 


P Si = 4 


_ a2 


xi 


giJ_ x 

6 + 2 


-4tA* 


+4ri 


L L 2 p z l 

y C os2P;L +(^ L -^-)sin2p J L 


■^■sin2p z L + 


i i 2 ft 


4ft 2 


cos2P/L- 


1 


4pi 


-2L 


P iL : 


sinPjLsinhPjL 


+4 u D xi 


-2L 


P iL : 


ft 


+b| 


sinpjLcoshPjL 

p^ 3 _i 

6 2 


(cos p / L sinh p z L ) + (sin P ,-L cosh p t L ) 

(cospjLsinhP/LHsinpjLcoshp/L) 
(cospjLcoshp J L)+(sinP J Lsinhp i L) 
(cosPjLcoshPjL)-(sinPjLsinhPjL) 


P, 


/ Z/ 2 p; i 

— cos2pjL+( - - 4^p~ )sin2Pi^ 


B x i Cxi 

+2L 


&L 4 


(sinP ,-L sinh P/ L )-(cos P, L cosh P, L ) 
cosP,LsinhP/L]--J- [ (cosP/LcoshP/L) 
+(sinpjLsinhPjL)' * 


-BxiDxi 

+2L 


P iL : 


P/ 


(sinp/Lcoshp,L)-(cospjLsinhpjL) 


cosPjLcoshp z L 


Pi 


(cospjLsinhpj-L) 


(24B) 
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P6i = 


+A 


+A 


-B 


-B 


+C 2 


XI 


2 


P,Lcosh2P,L+( 


-(sinPjLcoshp,L) 
1 


P,t 2 


2r 3 1 


4 P i 


)sinh 2 p,L- 


P i l L 


+CriD 


XL^XL 


+D 


XI 


•^~~cos2P,X-“sinh2P|L+'^~cosh2p i -L-“j- 


1 

2 


ft L 2 i R?L 3 ^ 

P,Lcos2p,L +( ~y~ + — )sinh 2 p ,L + ~y- 


1 yi 


P }L 3 j_ 

6 + 2 


AyiByi 


yiCyi 


L L 2 P i 

“C0s2p,L +( 2~ _ 4p~ )sin2Pi^ 


L 1 L 2 P/ 1 

fsin 2 M + (— )cos2 Pi L-— 


(25) 


Pit 1 


(cosp,Lsinhp,L)+(sinPjLcoshP J L) 


Pi 




P iL : 


-2 L 


sinP,Lsinhp,L 


(cosP jL sinhp,L )-(sinPjLcoshp,L ) 


(cosPjLcoshp,L)+(sinP J LsinhP,L) 


_1_ 

Pi 


-2L 


(cosPjLcoshPj-L^sinp/LsinhPjL) 


sinP,Lcoshp,L 
J_ 

Pi 


+5 




nCyi 


&L‘ 


6 2 


L L 2 3 ; 1 

— cos2p { L+(~^— )sin2p I L 


(sinP,Lsinhp,L)-(cosP,Lcoshp,L) 


+2L 


Pi 


(cosP,Lcoshp,L)+(sinP,Lsinhp,L) 




Pi^ 


(sinP,Lcoshp,L)-(cosP,LsinhP,L) 


+2L 


cosP,Lsinhp,L 

Pi 

cosP,Lcoshp,L 


Pi 


(cospjLsinhp/L^sinP/LcoshPjL) 


C 2 

'■'yi 


1 

2 


3 1 


pjLcosh2P,L+( 


fct" i ... . 2 „. Pit 


4 Pi 


)sinh z p,L-‘ 


CyiDyi 


B,L 2 t i 1 

cos2pjL-Ysinh2p,L+— COS2P/L-— 
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1 

2 


P;LC0S2P;L+( 




)sinh 2 P/L + 


3 



(26) 


The equations (25) and (26) can alternatively be derived by replacing P,- = 


a*. 

L ' 


The kinetic energy of the reflector is 

T 2 = (1/2) m2ylV 0 ~m2yla(L)u+m2yld(L)-(\/2) m 2 ® T a T (L)a(L)tii 

+m 2 w T 5(L)d(L)+( 1/2) m 2 d T (L)d(L)H 1/2) ft T l 2 0. (27) 

where m 2 is the mass of the reflector and / 2 is the mass moment of inertia matrix 
of the reflector. The deflection vector d(L) at the mass center of the reflector is 
given as 


d(L) = 


u x (L) - r y Uy(L) 
u y (L) + r x Uy 
u x \L)r x + u y '(L)r y 


(28) 


and the position vector from the point of attachment to the reflector mass center is 
given by 


a(L) = 



Thus, 


(29) 


d(L) = 


ii x (L)-r y Uy(L ) 
u y (L) + r x iiy(L) 
u x \L)r x + u y \L)r y 


(30) 


The angular velocity of the reflector in the inertial co-ordinate system ft can be 
shown to be 


ft = (0 + 


u x ' 
• * 
Uy 

My 


(31) 
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The equation (27) can be simplified as 

7*2 = (1/2) m 2 ylVo- m 2 yla(L)to+ m 2 Vld(L)+(l /2 ) m 2 L 2 J to? +co 2 j 
+m2® T a(L)d(L)+(V2)m 2 X X MLWxj(L )<Mj+ 

i=i/=i 

£ £ fymyjVtilij 

/= U=1 


where 


+(l/2)P r / 2 P+(l/2)(o 7 7 2 co (32) 


= . (33) 

,i=l i=l j=l 

The kinetic energy of the shuttle, T 0 , is given as 

7 0 = (l/2)miV r V+(l/2) co r J/i]to (34) 

where m i is the mass of the shuttle and / 1 is the mass moment of inertia matrix of 
the shuttle. 

The total kinetic energy is given as 

T = T 0 +T l+ T 2 (35) 


This can be simplified as 


T = (1/2) m 0 V T V+w T | //]k+( 1/2) w r J / 0 ] co+pL £ < 7 ?+Y r a 

+ 0 ) r rd+ 03 r p+/n 2 V r d(L)+m 2 co r rd(L)+ 


m 2 ©'a(L)rf(L)+(l/2)m 2 X 


+(1/2) pj 1 2 P+(l/4) p XP5i<?? + 


where 
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1 0 -1 ! +(1/3) pL 


m 0 = m\+pL+m 2 

H = pL+m2 

'10 0 
0 1 0 


0 0 0 


r+m2a(L)+pLc 


+I 2 +J 2 ~pLrr~pLrc -m 2 rr-m 2 ra (L ) 


The term J 2 in this equation can be shown to be: 


/ 2 = «2 


(ry+L 2 ) -r x r 


r x L 


y 

~r x r y ( rl+L 2 ) r y L 
r x L r y L (r 2 +r 2 ) 


The total kinetic energy expression can be further simplified as 

T = (1/2) m 0 Vj'VW \ tfl V+(l/2) to 7 ’ I " /„! [ A , 1 ? 

+W 7 '[a 2 ]£+( 


where 


= a +m 2 d(L) 

A 2J ^ = ra+$+m2rd(L)+m2a(L)d(L ) 


H = 


0 

pL+m 2 +P5i+P6i 

0 


In this equation 


(37) 


[ [<0 ( Z .)] r / 2 (<K ( i )]] . 


< mr = 


<t>i«U) 0 0 

0 *,/(£) 0 

0 0 4>, ¥ (L) 

teU) 6 6 

0 fyjyXL) 0 

0 0 fc ¥ (L) 
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Here i=2,3, ,n. The number n indicates the total number of flexible modes con- 

sidered. 


Equations of motion 

Lagrange’s equations of motion for the case of independent generalized co- 
ordinates qic are 


d dT dT „ dU 

— = — 

dt dq k oqic oq k 


(k = l,2, ,n) 


(38) 


where, T = T(q,q) is the kinetic energy 
U = U (q) is the potential energy, and 

Q k are the generalized forces arising from nonconservative sources. 


The generalized co-ordinates are: 

Rx,Ry> r z ~ position of orbiter mass center relative to inertial frame origin. 

01 , 02,03 - roll, pitch and yaw angles of orbiter. 

q\,q2, An ~ modal deformation co-ordinates for the beam. 


The previous kinetic energy expression developed in equation (37) is given 
in terms of nonholonomic velocities V and co, and generalized velocities q. Using 
the notation TCU.m,^) for this kinetic energy expression and T for kinetic energy 
expression in terms of generalized velocities, the equations of motion are 
developed. Thus, equation (37) is rewritten as 


T = (1/2) m„V T V + 0 ) r [tf]v + (l/2)co r [/„]o) + V T (37) 

+ 0) 7 [/l2]? + (l/2)£ r [A 3 ]4 


(a) Translational Equations 

From the chain rule applied to equation (37) using equation (8), one gets 
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dT 


dT 

dRx 


avi 

dT 

= C T 

ar 

dRy 


dV 2 

dT 


dT 

N 


dVi 


Also, the generalized forces are CF_{t) where 


( 39 ) 


Fit) = F 0 (t)+F 2 (t ) 


(40) 


Fp(t) represents the force applied at the orbiter mass center and F 2 (t) represents 
the force applied at the reflector mass center. From Lagrange’s equations 


_d_ 

dT 

^ • r 

ar 

dt 

dv 

+ CC 

dv 


= F(t) 


(41) 


and from equation (37) 


ar 

dV 


= m 0 V-H(£>+A\q 


(42) 


Substituting equation (42) in (41), 

m 0 V-H(n+A i'q = -CC T (m 0 V-H(0+A x q) + F(t) 


(43) 


This can be rewritten as 


m 0 V-Hw + Ai'q =Ni +F(t) (44) 

where the nonlinear term N i is given as 

N ! = -CC T (m 0 V-Hco+A ! q) (45) 

= -jb(m o V-H( 0 +A i q) 

Here, w - CC T . 

(b) Rotational Equations : 
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From equation (4) 


Again using the chain rule 


Also 


It can be shown that 


dV T 

30 f 


(0 


= M r 0 


dT 

30 


= M 


3 T 
3co 


3r 


‘ dV T 



3(o r 

30i 


30i 



30i 

3 T 


dV T 


dr 

1 

3(o T 

302 


302 


dv 

-r 

302 

3r 


dv T 

v — 


3co r 

303 


ae 3 



303 


= V'C 


t^SC t 


30,- 


.i =1,2,3,. 


and 


3(O r r.,-i 3 M . , n „ 

aST"^" mT ,=1 - 2 ' 3 ’- 


and 


3 T 
3co 


(46) 


(47) 


(48A) 


(48B) 


3 T 
30 


7>3C r 


- 30! 

dC T 


V l C 
V T C 
V T C 


- 302 

3 C T 
303 


31 

3V 


Tw-i 3M 
(o'M A — — 

- 30i 

t w-i 3M 
— 

302 

Tw-1 3M 
(0 y M * — — 

— 30 3 


3T 

3(0 


( 49 ) 
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From equation (37), 


and as before 


07 

3co 


= HV + 7 0 co + A 2 q 


(50) 


07 

dV 


m 0 V -He o + A\q 


Using the Lagrange’s equations 


(42) 


d_ 57 

* ae 



(51) 


where G is the net moment about the mass center of the orbiter with respect to the 
body-fixed frame. It is given as 


G = G 0 + ( r+a)xF 2 (52) 

Gp is the external moment applied about the mass center. Eqation (51) can be 
simpified by substituting equations (42), (49), and (50) together with the relation- 
ship developed in (46) as 


HV + I 0 w + A 2 q = G+ N 2 
where the nonlinear term N 2 is given as 

0C r 


— AA~ 1 


N 2 = M 


- ae, 

VTC ^~ 

— 802 

v t c SC t 


— 803 


07 

dv 


+M 


-1 


Taa - i 5M- 


co'M 


001 

Tt/i-llM- 


orM 


009 


~M\ 


T.f - 1 0M 

orM 1 


003 


(53) 


07 

00 ) 


. (54) 


(c) Vibration Equations of the Beam 
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Since T in equation (37) is given in terms of q which is a vector of general- 
ized velocities, 

37 _ 37 
3 q 3 q 

and 


3 7 

• 

The potential energy in the beam 


= AjV + + A$q 

is given by 


(55) 


U = (l/2)q T Kq 

where the stiffness matrix K is given as 


and 


K = 


\ 0 
ka 

0\ 


(56) 


(57) 


kn = EW 


J<j>i(5)ds + l^(j)d5 
0 0 


Li 

6* y Pyj J $ yj ( s )ds 


Gy represents the modulus of rigidity of the beam and Py,- = 

D is the mass per unit volume (mass density) of the beam. Thus, 


Dcof 


r y 


where 


3 U_ 

3 q 




( 58 ) 


Using the Lagrangian Equations (38) and assuming that F 2=0, 
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A\V + A 2 w + A^q = -Kq , (59) 

(d) Motion Stiffness 

As noted in references (27,28), generalized active forces associated with 
motion stiffness must be taken into account in order to obtain dynamical equa- 
tions for performing large arbitrary maneuvers. This provides the mechanism of 
compensating errors caused by premature linearization during the modal analysis 
for an arbitrary flexible body undergoing the most general large rotation and 
translation. 

For a mass element on the beam from the point of attachment (origin of the 
body-fixed frame) before deformation, the translational velocity in the inertial 
frame is as given as 


V a = V Q +(OX(a+d) + d(z,t ) (60) 

where V 0 is given by equation (9). The corresponding angular velocity in the iner- 
tial frame of the mass element is given as 

m a =03 + P(z,r) - (61) 

where 



0 

0 

Uy(z,t ) 


and «y(z, 0 is the torsional displacement of the beam. 

Differentiation of the right-hand side of equation (61) with respect to t yields the 
acceleration of mass element in the inertial frame as 


a a = V 0 + ri)x(a+d) + d(z,t) + to x V 0 + (ox(a+d)+2d(z,t) 


(62) 
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Similarly the angular acceleration of the mass center in the inertial frame follows 
from equation (4). 

a a = 0 ) + co x p(z,t) + p(z,t) (63) 

Equation (62) can be simplified as 

a a = V + cox(co xr) + wx(a+rf) + d(z,t ) 

+ cox [ V + coxr + co x(a + d) + 2<i(z,/)j (64) 

The inertia forces associated with acceleration in equation (64) can be given as 


F* = -J(p dz)a a - m 2 a a {L) 

0 

and the inertia moment associated with angular acceleration is 


(65) 


L L 

G* =-S(dI)a a -I 2 a a (L) -I 
0 0 


* 

(O a x [ W co«] -c o a (L)x j^/ 2 co] (66) 


The inertia force obtained in equation (65) is added to the right hand side of equa- 
tion (44) and the inertia moment of equation (66) is similarly added to equation 
(53). 

(e) Slewing Equations 

If it is considered to perform a slew maneuver about an arbitrary axis y and 
the slew angle to be then the slew maneuver can be expressed in terms of four 
Euler parameters. These four Euler parameters are defined as 


e = 


ei 

£2 


= X sin^- 

2 


(67) 


£4 = COS' 


( 68 ) 
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and their derivatives with respect to time are given as 


d£ 

dt 


= ( £4(0 + exco) 


d£ 4 
~dt 



( 69 ) 

( 70 ) 


The four Euler parameters can be related to the angular velocity components 
of the rigid assembly as 


■ -1 
El 


£i £4 -£3 £ 2 ‘ 


0 

k 


£2 £3 £4 -£1 


COi 

k 


£3 -£2 £1 e 4 


0)2 

£4 


£4 — £1 -£2 -£3 


(o 3 


( 71 ) 


If a slew maneuver is considered to be purely rotational, then the transla- 
tional velocity and acceleration can be shown to be negligible during the slew 
maneuver and only the rotational and vibration equations are reqired for the 
analysis and they are simplified by setting V= 0 in both ( 53 ) and ( 59 ) and are writ- 
ten as follows 


I 0 & + A 2 q = G(?)+N 2 ( co) (72) 

Ala + Aiq+Kq = Q(t) . (73) 

where, 

G(t) is the net moment applied about the mass center of the orbiter and is 
given by the following equations (figs. 3 & 4 ) 

G(t) = G 0 (t) + (r + a)xF 2 . ( 74 ) 

Also, £ 2(0 represents the generalized force vector which is given by the following 
equation 
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m 


Z (Qjx,(0 + Qjy\ (0) + 0^ + + 0^, 

7=1 
m 

£ (Qjx 2 ( 0 + Qjy 2 (f)) + 0* 2 + 0y 2 + Qyi 

7=1 


0(0 = 


m 


X ( Qjxi(t) + QjyX 0) + 0X1 + 0yi + 0 V* 

7=1 


(75) 


where, the generalized force components are given as 


L 

Qj Xi = l^(V)5(z-Zj)<M0<k (76) 

o 

L 

Qjy t = \f jy (z,t)§(z-z j)§ yi (z )dz (77) 

0 


and 


0, Vi (O = 0 . (78) 


Here, F jx (z ,0 is the x component of the concentrated force applied at location j 
on the flexible antenna and F iy is the y component of that force. 


Also, 


Qxt( I) = 

Qyid) = (79) 

0 y 0 ) — M y(f )§yi(L) 

Here, F 2 is the force applied at the reflector C. G. 

Thus, 
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M v (r) = F 2 x r y + F 2 y r x + M 2 V (80) 

The location of reflector C. G. is given by coordinates ( r x ,r y ) and M 2 ^ 
represents the external moment applied at the reflector C. G. 

Thus equations (79) - (80) completely represent the dynamics of the slew 
maneuver. These equations are nonlinear and coupled including both the rigid- 
body dynamics and the dynamics of the flexible appendage with kinematic non- 
linearities. It is important to note that the nonlinear term A^ 2 (co) is dependent on 
the rotational velocity and as a result determined by the slew maneuver rate. Thus 
the basic slew maneuver stretegy has to be developed before this term can be 
linearized. 

(f) Vibration Equations of the Beam with Damping 

If damping is included in the derivation of vibration equations of the beam, 
then the damping effect can be expressed in terms of frictional forces. These are 
nonconservative, retarding forces and are assumed to be proportional to the gen- 
eralized velocities. In deriving the vibration equations by means of Lagrange’s 
equations, the following function is introduced 


Fd = T £ £ bijWj . (81) 

z i= v-i 

It also has a positive definite quadratic form similar to the kinetic and poten- 
tial energy expressions. 


With this definition, Lagrange’s equations assume the form 


d dT 

dtdq k 


dT d?d 
dqk 


= Qk~ 


dU 

dqk 


(£=1,2, ,n) 


(82) 
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Again, as before 


and 


dT 

dq 


A\v + aJco + A 2 q 


(55) 


dU 

dQ 



(58) 


and it can be seen from (81) that 




where the damping matrix B is symmetrical and is given as 


(83) 


'*>11 

b\2 • 

• *>l/i 

*>21 

*>22 • 

• *>2n 

b n \ 

*>n2 ■ 

• *>nn 


(84) 


The vibration equations are given as 


A\V+Al® + A 3 q+Bq=-Kq + Q(t) . (85) 

The slewing equations (72) and (73) would be modified as 


/ 2 d) + A 2 q = G(0 + N 2 ((0,q) (86) 

Al& + A 3 q+Bq =-Kq + Q(t) . (87) 

Nonlinear Term in the Rotational Equations 

The nonlinear term N 2 in the rotational equations (72) and (86) during the 
slewing maneuver is simplified as 
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where 


( 88 ) 


r > ,-i OM 

V" 


- = [ooo 


(89) 


r, f _i OM 1 

or M 1 — - = — 

— 002 COS02 


^ (-(Oisin02cos 2 03+co2sin02sin03cos03) (coisin02sin03cos03 


-(O2sin02sin 2 03) ((Dicos02cos03-ot)2cos02sin03) 


(90) 


f ((O2COS02) (-COiCOS0 2 ) 0 ] (91) 

— 003 COS02 L J 

Since the transformation matrix, M , is a function of 02 and 03, the time 
derivative of M can be expressed by the chain rule as 


OM \ dM i 

M = + 3eT 03 


(92) 


From equation (5) 



(~sin02cos03)02 (sin02sin03)02 (cos02)02 
0 0 0 
0 0 0 


(93) 


OM 

003 



(-COS 02 Sin 03)03 (-COS 02 COS 03)03 O' 

(COS03)03 (-sin03)03 0 

0 0 0 


(94) 


Substituting these equations (93) and (94) in (92) 
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M = 


(-sin 02 COS 03 ) 02 +( - CO^ 02 s in 03)03 (sin 02 Sin 03 ) 02 +(-COS 02 COS 03>03 (COS02>02 
(cos03)03 (-sin03)03 0 

0 0 0 

(95) 


From equation (4), this can also be expressed as 


M = 


COS02 


(-sin02cos03)(coicos02sin03 (sin02sin03)(coicos02sin03 

+CO2COS02COS03)+(-COS02Sin03) +CO2COS02COS03)+(-COS02COS03) 
(-a>isin02cos03+a>2sin02sin03 (-(Oisin02cos03+co2sin02sin03 

+CO3COS02) -KO3COS02) 

(co^ 03 )(-a)isin 02 cos 03 (-sin 03 )(-a)isin 02 cos 03 

+(O2sin02sin03+co3COS02) +co2sin02sin03+a)3cos02) 


0 


0 


COS02((OiCOS02Sin03+CO2COS02COS03) 


0 

0 


(96) 


Also, M 1 is given as 


M~ l = 


1 

COS02 


cos©3 cos02sin03 -sin02cos03 
-sin03 cos02cos03 sin02sin03 
0 0 COS02 


Thus, the nonlinear term N 2 can be rewritten as 


( 97 ) 


Where the term A '3 is 


N 2 =A' 3 (co,0) I 0 <o+A 2 q 
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A' 3 (co,0) = M 


- ii/f-i 


0 

t. ,-i 0M 
(Sr M 1 


ae 2 


t. i 0 m 
co'M 1 


ae 3 


-Ml 


#2 = A' 3 (co,0)/ O co + A' 3 (cd,0)A 2 <? 


= A 4 (<n,0) + A 5 (co,0)<7 (98) 

where A 4 depends on the rigid-body slewing and is nonlinear in terms of co and 0. 
The second term relates the coupling between the rigid-body slewing and the 
flexible modes. This equation can be further simplified in terms of Euler parame- 
ters by relationships developed in appendix II as 


N 2 = A 6 (co,e) + A 7 (co,e)<? (99) 

where £ is the Euler vector comprising all four Euler parameters, 

From equations (72) and (73) and by defining A = A 2 T 1 0 ~ l A 2 + A 3 , the 
following equations are obtained 


(b = / 0 1 


A 2 A ^ Bq +A 2 A ^ Kq + 


-l, 


A 2 A 1 A 2 t I 0 1 +h 


G(t) 


1 

■ 

a- 1 'a 1 t i 0 -' + / 3 

W 2 (c o,e) 


q = A~ l Bq - A~ l Kq - A-'A 2 T I 0 ~ l G(t) + A- l A 2 T l 0 - l N 2 ( co,e) 


( 100 ) 


+ A- 1 G(r) 


( 101 ) 


It is assumed that control forces applied for vibration suppression has negli- 
gible effect on rotational maneuver of the spacecraft in developing equations 
(100) and (101). Also, / 3 represents 3x3 identity matrix in these equations. 
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Subsystems and State Variable Models 

The two dynamical subsystems considered for decentralized control are the 
dynamics of the slewing of the rigidized SCOLE assembly and the vibration 
dynamics of the flexible antenna. These subsystems are represented by subscripts 
I and II respectively for subsequent analysis. 

The following are the definitions of state variables and control variables for 
subsystem I. 

xhAei\ x /2 A £2; x 7 3 A £3 ; x/ 4 Ae 4 ; 

xjs A ©1; */6 A CO2; Xj 7 Aco 3 ; 
x/gAGi; X79AG2; X/10AG3; 

«/i AGi; uj 2 A G2> m/3 A G 3 . 


The interconnecting variables from the second subsystem to reflect the cou- 
pling between the subsystems are defined as 

z/iAr//i; Z12AX112', Z/3AX//3; z /4 Axjj 4 \ zj 5 AO; 

2/6 AO ; Z/7AO; z/8 A 0 ; z /9 AO ; z /10 AO. 


The following state equations are obtained for subsystem I using these 
definitions 


xi 1 
Xl 2 
xn 

f/4_ 

Xl5 

Xl6 

_XlJ__ 

XJ9 

*/10 



*/l ' 



xi 2 



X I3 


o| 0 j_0 

Oj_Z); 0 


_X_I4_ 


o' l o[h_ 

1 1 


XI 5 
X I6 

+ 

1 1 

0 1 Z?2 1 ^3 
, - - 

0» 0 « 0 


Jill _ 


0 

0 


""•*78 



Xl 9 



*/l0 



M/1 
«/ 2 
M/3 
2 /l' 
.2/2 
2/3 
2/4 
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B 2 =Io~ X A 2 A- l K , 

b 3 = i 0 - 1 a 2 a~ 1 b , 

Xi 4 -Xi 3 x/ 2 
Xi 3 X/ 4 — X/ 1 
-X 72 X/i X /4 
“*71 -X/2 ~X] 3 

For the second subsystem which is the flexible appendage of the entire sys- 
tem, the first two flexible modes are considered and the corresponding state vari- 
ables and control variables are defined as 




Xn\Aq\', Xjj 2 A q 2 \ 
xmAq\\ Xu\Aq 2 \ 
uinAQi; u II 2 AQ 2 \ • 


As in the previous case, the coupling between the two subsystems is derived 
in terms of the following interconnecting variables from the first subsystem. 

Z 771 A X/ 1 ; Z //2 A X/ 2 ; Z 7/3 Ax/ 3 ; 

Z//4 Ax/4; 


Z 775 AX/ 5 ; z // 6 Ax /6 ; Z// 7 AX/ 7 ; 
z /77 A x/g; z //9 A x/ 9 ; zjj 10 A x / 10 . 
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The following are the state equations of this subsystem. 


• 



Xjl 1 


_ l 


Xu 1 

xin 


0 ■ I 


Xll 2 

XII3 


-A -1 K ' -A" 1 B 


~XU~3 

Xll A 


1 


Xll A 


o 

0 


' -a- 1 a 2 V‘ 



M//1 


-HI 


2 11 8 


2 II 9 


. z //10. 


+ 



-A-U/V 1 


N 2 ( z // 1 » z /7 2 > z /7 3 » z /7 4 > z 77 5 > z 77 6 > 2 II 7 ) 


The Optimal Control Problem 


(103) 


A general problem for the optimal control of interconnected dynamical sys- 
tems like large flexible spacecrafts can be formulated as 


Minimize J (Xj,Wj,z,) i = \,2,...,N (104) 

w.r.t. uj 

where x,- is the n, dimensional state vector of the ith subsystem, «,• is the 
corresponding m t - dimensional control vector and z, is the r, dimensional vector of 
interconnection inputs from the other subsystem. The integer N represents the 
total number of subsystems and the scalar functional J is defined by 


N 

i=l 


Pi ) + j Li [ x,(r),u,(r),z,(o] dt 

to 


(105) 


where Li |^x,(f),M,(0» z /(0j is the performance index at time t for i = 1,2,..,N sub- 
systems. The functional J defined in equation (105) is to be minimized subject to 
the constraints which define the subsystem dynamics, i.e. 


=fi Xi(t),Ui(t)Mt),t ,t 0 <t<t f 
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Xi(t 0 ) = x io , i = 1,2 ...,N . 

Also, the minimum of J must satisfy the interconnection relationship 


( 106 ) 


£ G* (Xi,Zj ) = 0 . (107) 

j=i 

The Open-loop Hierarchical Control 

Using the method of Goal Coordination or infeasible method [17,22], we 
consider another problem which is obtained by maximizing the dual function 
A, ) with respect to X(t) ( t 0 < t < tf ), where 


<P( X, ) = Min 
x,u,z 


J(x,u,z,X) 


(108) 


subject to constraints in equations (106) and (107). Here 


X = 

X\ 

u = 

Ml 

Z = 

£l 


X N 


U_N 




(109) 


Also, A, in equation (108) is a vector of Lagrange multipliers which is given as 


X = 


h 




( 110 ) 


N 
i = l 


Pi ( Xi(tf ) ) + | Li ( Xi,Ui,Zi,t ) dt + ]x T c i* ( Xi,Zj,t ) 


dt 


( 111 ) 
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j = U2,...,N 


Rewriting this functional J as 


N 

2 

1=1 


N 

J = T,Ji 

i=l 


Pi( Xi( tf)) + j Li ( Xi,Ui,Zj,t ) + }JGi{ Xi,ZiJ )1 dt 

to 1 J 


( 112 ) 


where. 


Ji = Pi( *i( t f ) ) + j [ Li ( Xi,Ui,Zi,t ) + }jGi ( Xi,Zi,t) j 


dt 


(113) 


and where 7J G*( Xi,Zj,t ) has been refactored into the form )J Gi(Xi,Zi,t ), i.e. 
into a form separable in the index i. 


Thus, 


J =J + 


t 0 ‘= 1 


(114) 


Then by the fundamental theorem of strong Lagrange duality [25] 


min J = max O ( X ) , i = 1,2,...,1V. (115) 

Uf X 

Thus an alternative way of optimizing J is to maximize O ( X ) . 

From equation (112), for a given X ( t ) , t 0 < t < tf , the functional J is separ- 
able into N independent minimization problems, the ith of which is given by 


Min Ji = Pi ( xi (!/)) + / Li ( Xi,Ui,Zi ) + XFGi ( Xj,Zj ) 1 

t J 


dt 


(116) 
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X h Ui,Zi 

subject to 



This leads to a two-level optimization structure where on the first level, for 
given X y the N independent minimization problems described in equations (116) 
and (117) are solved and on the second level, the X (t) ( t 0 < t < tf ) trajectory is 
improved by an optimization scheme like the steepest ascent method, i.e. from 
iteration j toy +1 


X 0 fV +1 = X (ty + a) + d j ( t ) t Q <t<tf (118) 

where 


di =V<*>(X(t))= XGiiXiyZj) , (119) 

i=i 

V ( X) is the gradient of <J> ( X) , a 7 > 0 is the step length and dJ is the steepest 
ascent search direction. At the optimum dJ — > 0 and the appropriate Lagrange 
multipier, X , is the optimum one. 

The development of this algorithm depends on the assertion Max d>( A.) = 
min / and this may not be valid for all nonlinear systems. Consequently, lineari- 
zation of Gi , and linearized equations for^- may be required for constraints to be 
convex and convexity of the constraints is necessary to prove this assertion. 
Nevertheless, the method is attractive from the standpoint of simplicity and that 
the dual function is concave for this nonlinear case. This ensures that if the dual- 
ity assertion is valid, the optimum obtained is the Global Optimum. 
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On the first level, since equation (116) is to be minimized subject to equation 
(117), the necessary conditions lead to a two point boundary value problem from 
which an open loop optimum control could be calculated. However, it is desirable 
to calculate a closed loop control and for this the quasilinearization approach can 
be utilized at level one for all subsystems. Thus an iterative scheme can be set up 
whereby an initial trajectory of X (t)* , t 0 < t < tf is guessed at level two and pro- 
vided to level one. At level one the two-point boundary value problems of the 
subsystems are solved by quasilinearization. The state and control trajectories of 
all the subsystems obtained at level one are sent to level two. The test for 
optimality based on equation (119) is conducted at level two and if this is not 
satisfied, equation (118) is used to obtain the new X (t) for the next iteration. 

Subsystem Closed Loop Controllers 

The closed loop controllers are obtained at the first level by solving the two- 
point boundary value problems of the subsystems utilizing the quasilinearization 
procedure. As noted in equations (116) and (117), the first level problem for the 
ith subsystem is 

For given X(t) , t 0 < t < tf , 


mm 


Xi,U h Zi 


p i [ Xf( <y)] + I [a( Xi,Ui,Zj ) + X T Gi( x h zj) j dt 


subject to 


Xi ~fj (Xj, U h Zi ), t 0 <t<t f 


(116) 


Xi (to) — Xio • 


( 117 ) 
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For this problem, the Hamiltonian //,• can be written as 


Hi = Li ( Xi,Ui,Zi ) + X r Gi ( x it Zj ) + r\ i T fi ( Xi,Ui,Zj ) . 
For a given X, the state and costate equations become 


Xi(t)=fi(Xi,Ui,Zi) 


r\i (O = - 


dHj 

dxj 






with 


( 120 ) 


( 121 ) 

( 122 ) 


dHj 

dui 




(123) 


It is assumed here that using the equations (122) and (123), it is possible to 
obtain the control «,• and the interconnect variable vector z, which is an explicit 
function of Tj i and jq, i.e. 


Hi = £i (*/,% ) 

(124) 

^i = di ( ) 

Using these relationships for «,• and z { in equations (121) and (122), the following 
equations are obtained 


Xi =Oi (x,-,T|j ) , 

t 0 <t<t f 

(125) 

r\i =bi (£/,%), 

VI 

VI 

(126) 


with the boundary conditions 


%i ( t 0 ) — Xjo 

and from the transversality conditions 


(127) 



(128) 


n /(//) = 


dPi 


Xi ( tf ) 


dxi 


Quasilinearization Procedure 


The two-point boundary value problem of ith subsystem is given by equa- 
tions (125) and (126) subject to boundary conditions of equations (127) and (128). 
This problem is solved by quasilinearization technique as follows. 

Defining y = 

equations (125) and (126) can be rewritten as 

y(/) = f[>(()] ■ (129) 

In the quasilinearization procedure, starting from an initial guessed trajec- 
tory for y = yi (t) , successive linearizations are performed of equation (129) in 
such a way that the final linear equation fory solves equation (129) to an accept- 
able degree subject to boundary conditions (127) and (128) which could be 
expressed in a more general form as 


y ( ‘ 0 ) t a 0 = bj 

(130) 

y ( t f ) T A f = b/ 

(131) 


where A 0 , Af are 2n x n matrices. 

The linearized equation of (129) about a trajectory y = y 7 (t ) is obtained by 
Taylor series expansion as 



y =F(y j ) + J(y j )(y-y j ) + 'v 


(132) 


where J (y J ) is the Jacobian of F 


y(t) , t 0 <t <tf, at y J and ¥ represents the 
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contribution of the higher order terms. Neglecting these higher order terms, the 
following linear equation is obtained 


y_ = F(yJ)+J(yj)(y-yj) . (133) 

If the initial guessed trajectory y-> while satisfying equations (130), (131) and 
(133) does not satisfy equation (129), then an iterative search can be utilized to 
obtain a better linearizing trajectory by various methods discussed in references 
19, 20, and 21. This iterative search is given by noting that equation (133) can be 
written by expanding individual equations (118) and (119) by Taylor series 
expansion about a known trajectory xj(t), v\{(t), t e [ t 0 ,tf] , and retaining terms 
of up to first order. The linearized reduced differential equations are 


it 0+1) = q_i [ V (0. V (0 


+ 


3 ( 2 ; 

~T~ (x{ ( t ) ,r\{ (0) 

dXi ~ — 


3a,- 

^i(0) 

3rj, - 


x,°' +1) (0-x/(0j 

^ +,, «)-il/«)] (134) 



Xi j (t),T]i j ( t ) 

l 4 

+ 


+ 

dbj 

(xj ( 


( x{ (o ,ni (0 ) (o - x{ (/)j 

UP*” «)-])/(()] (135) 


These differential equations can be rewritten as 1 


Xf m) =^u«)s 0 ' +1) (0 +^i 2 (Oni 0+1) «) + £i / (0 (136) 

lii 0+I) = A 2i <.ty S V+'\t) + A 22 (I)T| i 6 +1 >(r) + ejif) (137) 


or, in the partitioned matrix form, 


' ii 0>1) (0' 


ii/ 0+1> a) 



^ ll (0 ^ 12(0 


’ XjV +l Xt) 


£ i J « 

^ 2 l (0 ^ 22(0 


ri ,° +1) (0 

+ 

ei'(t) 


( 138 ) 
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where the matrices 


da; da; 

A U(0 A-, A 12 (OA — , 


db{ dbi 

A 21 «)A — , a 22( ,)4^7. 


and 


£i j A - A 11 (t)Xi j (t) - A nitWit) + a t 

(139) 

e 2 J a - A 2 i(t)Xi J (t) -A 2 2(tW(t) + bi 
are evaluated at Xj J ( t ), r\i J (t) and hence are known functions of time. 

The method of complementary functions [24] can be incorporated with this 
linearization of the differential equations in the implementation of iterative 
search. 

An initial guess, xf, T|,°, t t[t 0 , tf], is used to evaluate matrices in equa- 
tions (139) at the beginning of the first iteration. In the next step, n sets of solu- 
tions to the 2n homogeneous differential equations 


ii 0+1) =Ai 1 (<fe 0 ' +I) «) + A 12 (/)ii i 0' +1 )(r) (140) 

% c,+1) = + AjjfijnMo (i4i> 

are generated by numerical integration. For (J+ 1) st iteration, these solutions are 

denoted by x/* 1 , x/ 12 , t| i H2 \ \x^ n , r\ , Wn . The boundary conditions 

used in generating these solutions are 


Xi H1 ( < o ) = 0. 

n i H \ t 0 )= [ 100 . 

..0 

Xf H2 (. to ) = o , 

tf\to)= [oio . 

. .0 
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(142) 

X/ H "( > 0 ) = 0 , %"”(<<,)= [ooo ...l] 7 . 

Next, one particular solution at 0 + 1) denoted by xf , r | t p , is generated by numer- 
ically integrating equation (138) from t 0 , to tf, using the boundary conditions 
xf(t 0 ) = x 0 , T\i P (t 0 ) = 0. Then, the complete solution of equation (138) can be 
obtained by using the principle of superposition and is of the form 

Xi^ +l \t) = C\Xi Hl (t) + C 2 Xj H2 (t) + .... + C n Xj Hn (t) + xf(t) (143) 

1 1i° +1) (0 = C + c 2 T\i H2 (t) + .... + c n r\i Hn {t) + T\i P (t) (144) 

where the values of c i , c 2 , . ... ,c n which make q/ 1+1 (t/) = r\f are to be deter- 
mined. To find these values of c i , c 2 , . . . . , c n , we let t = tf in equation (144) 
and write it as 

T l/= []1 i Hl (t f ) r\i H2 (t f ) . . . l c+r\i p (t f ). (145) 

Here, c A ^ c i c 2 . . . . c„j is unknown. Solving for c yields 

£ = [ T]i H 1 (t f ) i ]i H2 (t f ) .... %"%)] 1 [ 1 y “ ^ ■ ( 146 > 

It is important to note that the indicated matrix inversion in equation (146) 
has to exist in order to solve for c. Substituting this solution of c into equations 
(143) and (144) gives the O'+l) st trajectory. This completes one iteration of the 
quasilinearization algorithm and this trajectory can be further utilized to begin 
another iteration, if required. Generally, the iterative scheme is terminated by 
comparing the j th and j + 1 st trajectories by calculating the norm shown in the 
following equation and comparing it with a preselected termination constatnt, p. 



( 147 ) 


Closed Loop Control 


*.0 + D 

T| l ^ +1 ) 





In order to obtain the closed loop control, the solution of the linearized equa- 
tion (133) can be written as 


u 

y_{ */) = < H tf,t)y_ (t) + J <j> (tf, x) [ F (x) - J y (T) j dx (148) 

where <J) is the state transition matrix of the system in equation (133). Rewriting 
equation (148) in terms of solutions of states and costates and replacing the 
integral terms by p,- (t) 


Xi Of) 


<>n Of, 0 <>12 Of, 0 


Xi (0 

1 

Pil 0) 

11/ Of) 


<i>21 Of, 0 <l>22 Of, 0 


11/ (0 

T 

Pi 2 (0 


From equations (128) and (149) 


(149) 


% Of) = 


dP)_ 

dxi 


= <t>21 Of, t)Xj {tf, t) + ({>22 Of, 0*1/ (0 + Pi 2 (0 • 


(150) 


Thus, 


11 / (0 = 


dPi 

— 4>21 Xi (0 — P/2 (0 



dpf 

dxi 


~Pi 2 (0 




(151) 


It is important to note here that <j> 22 -1 always exists since it is a principal 
minor of the state transition matrix. 

Substituting equation (151) into equation (124) 
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Ui = Ci(Xi,t) . (152) 


4. NUMERICAL DATA 


The analytics developed in the previous section are utilized together with the 
basic SCOLE data [1] and the three dimensional linear vibration analysis [4] to 
generate the following numerical data. 


m i = 6366 A6slugs. ; m 2 = 12 Alslugs. ; p = 0.0955 slugs /ft. ; L = 130 'ft. 

= 72E+%lb/ft 2 \ C EI) X = ( EI) y = (El) = AE +llb-ft 2 \ 


' 0.036 ' 
r = -0.036 

“ -0.379 

' 0 
c = 0 

' . -65J 

905443.0 0.0 145393.0' 

I x = 0.0 6789100.0 0.0 

.145393.0 0.0 7086601.0 

' 18000.0 -7570.0 0.0 ' 

l 2 = -7570.0 9336.0 0.0 

0 0.0 27407.0 


The three dimensional vibration analysis is given in terms of the first ten 
modal frequencies and mode shapes in table 1 . Here, 

a ^ _ a is a.iS clis 

tyxi(s) = sin —— +Bxi cos — ~ sinh —— +D x ,cosh 

L L L Li 


a ^ a ^ a,s a is 

<{>^•( 5 ) = A yi sin —— -t-fiy,- cos —— + Cyi sinh +D yi cosh 


s s 

0\|/i(^) = A yj-sinOyi — +B y,cosa vi — 
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a-i = 


(ofpL 4 


1 

4 


El 


<V — 


DL 2 c o} 
G 


1 

2 


Using these data the following matrices are obtained. 



' 1216640 - 1.530307 175667.1 
- 31.66433 7082976 - 52474.84 
175690 - 52503.9 7131493 
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THREE DIMENSIONAL MODE SHAPE CHARACTERISTICS 


MODE No. 


FRE 


a 

A x 

B x 

C x 


MODE No. 


FRE 


a 

Ax 

B x 

Cx 


MODE No. 


FRE 


a 

Ax 

B x 

Cx 


1 

2 

0.27804240E+00 

0.31357296E+00 

O.12O12084E+O1 
0.16282665E+00 
-0.19670286E+00 
-0.16983450E+00 
0.19616259E+00 
-0.1 02746 18E-01 
0.57579 133E-02 
0.11810057E-01 
-0.57220462E-02 
0.19360955E-01 
-0.50748354E-01 
0.13978018E-04 

0.1 27565 18E+01 
0.38855291E-02 
-0.14998387E-01 
-0.43321018E-02 
0.14985820E-01 
0. 1421978 1E+00 
-0.22695797E+00 
-0.1 9283 105E+00 
0.226445 6 1E+00 
0.21835058E-01 
0.31115282E-01 
-0.75992337E-05 

3 

4 

0.81 300 189E+00 

0.11856099E+01 

0.20540387E+01 
0.40868 188E-01 
-0.61958845E-01 
-0.41309992E-01 
0.61880796E-01 
-0.22438404E-01 
0.36509234E-01 
0.24390447E-01 
-0.36464758E-01 
0.5661 1842E-01 
0.9269890 IE-01 
-0.87320799E-05 

0.24804687E+01 
0.80641794E-01 
-0.67233377E-01 
-0.80913938E-01 
0.671 063 16E-01 
0.13728679E+00 
-0.11746932E+00 
-0.14085209E+00 
0.11725057E+00 
0.82557693E-01 
-0.16158934E-03 
0.1 04377 18E-07 

5 

6 

0.20536300E+01 

0.497 16090E+01 

0.32645546E+01 
0.99278 129E-01 
-0.92344553E-01 
-0.99442 145E-01 
0.92225801E-01 
-0.57396019E-01 
0.5397 6008E-01 
0.581 14853E-01 
-0.53906980E-01 
0.14300062E+00 
-0.16588614E-02 
0.61861804E-07 

0.497 16090E+01 
0.45739784E-01 
-0.46365581E-01 
-0.45763 106E-01 
0.4632967 6E-01 
0.78612940E-01 
-0.79952853E-01 
-0.789 14485E-01 
0.79891039E-01 
0.33165303E+00 
-0.93394833E-05 
0.1501721 IE-09 
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THREE DIMENSIONAL MODE SHAPE 


MODE No. 


CHARACTERISTICS 


MODE No. 


0.55157833E+01 

0.12281249E+02 

0.53501560E+01 
0.81311804E-01 
-0.82056569E-01 
-0.81344923E-01 
0.81997259E-01 
-0.47145439E-01 
0.47703590E-01 
0.47289807E-01 
-0.47669 155E-01 
0.38408 110E+00 
-0.23855560E-02 
0.331 2204 IE-07 

0.79833305E+01 
0.44835061E-01 
-0.448349 14E-01 
-0.44840508E-01 
0.448 13000E-01 
0.77404756E-01 
-0.77465629E-01 
-0.77475327E-01 
0.77427782E-01 
0.8551 8 143E+00 
0.15830371E-05 
-0.987 15017E-11 

9 

10 

0.12890442E+02 

0.23679520E+02 

0.81789349E+01 
0.78743585E-01 
-0.78755259E-01 
-0.78752483E-01 
. 0.787 17693E-01 

i -0.45569244E-01 

0.45609474E-01 
0.45607884E-01 
-0.45587726E-01 
0.897 60145E+00 
0.94995483E-03 
-0.56437766E-08 

0.11 085 347E+02 
0.44348498E-01 
-0.44367373E-01 
-0.4435051 IE-01 
0.4435 1763E-01 
0.76707490E-01 
-0.76762782E-01 
-0.7673361 2E-01 
0.76735779E-01 
0.16488784E+01 
-0.51105957E-06 
0.16528495E-11 
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' 0.12927£+2 -0.46332£-l -0.56706£-1 -0.23671£-1 -0.24099£0 
-0.46332£-l 0.1293 IE +2 0.50635£-l -0.48234£-l 0.13228£0 

-0.56706£-l 0.50634£-l 0.13644£+2 -0.29552£-2 -0.30599£0 

-0.2367 1£-1 -0.48234£-l -0.29553£-2 0.12541£+2 0.44849£-2 

-0.24099£ 0 0.13228£0 -0.30599£0 0.448486£-2 0.1324£+2 

0.32322£-2 0.651 14£ -2 -0.41218£-3 -0.11644£-1 -0.22449£-3 

0.12271£0 —0.71 133£ — 1 -0.19374£-1 -0.14618£-2 -0.32102£0 
—0.1 1 158£ —2 -0.23852£ -2 0.22451£-3 0.42121£-2 -0.10488£-3 
-0.57067£-1 0.34303£ -1 0.65492£-l 0.40936£-3 0.11851£0 
0.12044£-2 0.10106£-2 0.54187£-3 -0.23616£-2 -0.18768£-2 


0.32323£-2 0.12271£0 -0.11158£-2 -0.57067£-l 0.12045£-2' 

0.651 14£ -2 —0.71 133£— 1 -0.23852£-2 0.34302£-l 0.10106£-2 
-0.41 21 8£ -3 -0.19374£-1 0.2245£-3 0.6549 1£-1 0.541 87£ -3 

-0.11644£-1 -0.14618£-2 0.42121£-2 0.40936£-3 -0.23616£-2 
-0.22449£-3 -0.32102£0 -0.10488£-3 0.11851£0 -0.18768£-2 
0.12432£ +2 0.17594£-3 -0.57146£-3 -0.10095£-3 0.31955£-3 
0.17594£-3 0.12606£ +2 0.82787£-5 -0.67559£-l 0.76582£-3 
-0.57146£-3 0.82787£-5 0.12425£+2 0.10429£-4 -0.11519£-3 
-0.10095£-3 -0.67559£-l 0.10423£-4 0.12477£+2 -0.2948£-3 
0.3l955£-3 0.765 82£ -3 -0.11519£-3 -0.2948£-3 0.12423£+2 



0.7142542£2 
0.2384121£ +3 
-0.5915277£ +2 
0.2843785£ +3 
-0.3872402£ +3 
-0.2010108£2 
0.1731334£ +3 
0.1495208£ +2 
-0.81759£ +2 
-0.3200308£ +1 


-0.3890922£ +3 
-0.3622244£ +2 
-0.1761796£ +3 
-0.469601 1£ + 1 
0.6004922£ +2 
0.3311526£2 
0.7672109£ +2 
-0.2067402£ +2 
-0.6396428£ +2 
0.1692286£ +2 


0.9536207£ + 1 ' 
0.27423 15£ + 1 
0.2077235£ 2 
0. 1 62303 9£ + l 
-0.1236509£ +1 
-0.1009266£ +1 
-0.47943 85£ +1 
0.642263 1£0 
0.3995 166£1 
-0.4793963£0 
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The stiffness matrix K is calculated using equation (57) and the mode shape 
coefficients given in Table 1. This matrix is a diagonal matrix and is represented 
in terms of the diagonal elements as 


k u = 0.3666744£2 ' 
k2 2 = 0.4647383£2 
* 3 ’ 3 = 0.3148351£3 
k 4 ’ A = 0.687065 IE 3 
k 5<5 = 0.2065265£4 
k 6t6 = 0.1114602£5 
k-j'j = 0.14901 14£ 5 
fcg’g = 0.7391932£5 
k 9>9 = 0.8131641£5 
k 10,10 = 0.2748996£6 


The damping matrix B used for this analysis is a diagonal matrix and for 
damping ratio £ = 0.003, it is calculated to be 


6 u = 0.9685964£-3 ' 
b22 = 0.1088608£-2 

& 3 ’ 3 = 0.28340 16£ -2 
b* A = 0.4256808£-2 
b 5 ] 5 = 0.7387177£-2 
be ,6 = 0. 1 7 1 90 1 4£ — 1 
b 7 'j = 0.1984237£-1 
Z> 8 ’ 8 = 0.4421234£-1 
b 9 \ 9 = 0.4633434£-l 
b mo = 0.8527647£-l 
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(c) 


(d) 



(a) Axes in reference position (b) First rotation-about x axis 
(c) Second rotation-about y axis (d) Final rotation-about z axis 


FIGURE 2 
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Figure 3- Position Vectors in Inertial Frame 
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Figure 4- Vectors in Body-fixed Frame 
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APPENDIX I 


The following is a summary of transformations between inertial frame and 
body-fixed frame. Here, q and q (/ =1,2,3) denote sinG,- and cosG; (/ =1,2,3) 
respectively. 


(a) Space-three Angles 


C = 


C2 C 3 C 2^3 ~S 2 

S\S2Cl~S2C\ S\S2S2+C2C i S \C 2 
ClJ , 2 c 3+S3‘ y l Cir2^3 _ C3‘y 1 C 1 C 2 


(b) Space-two Angles 


M t = 


1 0 -52 

0 Cl S\C2 
0 —5 1 C1C2 


C = 


C 2 S 2 S 3 -S 2 C 3 

S\S 2 S\C 2 S 3 +C 3 Ci S 1 C 2 C 3 +S 3 C 1 
C\Si -C1C2S3-C3S1 C1C2C3-53J1 


M 


T 


1 0 C 2 

0 c 1 5^2 
0 -si cis 2 


(c) Body-two Angles 


C = 


C 2 -*152 -C 1 S 2 

S2S3 —s 1C2S3+C3C 1 C1C253+C351 
S 2 C 3 -S 1 C 2 C 3 -S 3 C 1 C 1 C 2 C 3-5351 


M t = 


C 2 0 1 

5 2 ^3 C 3 0 

S 2 C3 S 3 0 
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APPENDIX II 


The transformation that relates the orientation angles 0 to Euler parameters 8 
is a nonlinear transformation. This transformation is developed for body-three 
angles representation in this appendix and similar transformations can be derived 
for other three representations, namely space-three angles, space-two angles, and 
body-two angles. 

(a) For sin 02 1 : 


If--?- < 02 < “.then 


02 = sin 1 


2(8381+8284) . 


If ( cos 0 icos 02 ) > 0 , then 


(A.l) 


01 = sin 1 


If ( cos 0 icos 02 > < 0 , then 


-2 (e 2 e 3 - 8184 ) 





cos 

sin 1 2( 8381 + 8284 ) 



01 = it- sin 


If ( cos 0 2 cos 0 3 ) > 0 , then 


-1 


-2 (8283-8184 ) 


sin 1 2( 8381 + 8284 ) 


cos 


(A. 2 ) 


(A. 3 ) 


03 = sin 1 


cos 


-2 (ei8 2 ~ 8384 ) 
sin" 1 J 2( 8381 + 8382 ) 


(A. 4 ) 
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If ( COS02COS03) < 0, then 

-2 (£ i &2 - £ 3^4 ) 
cos sin -1 ^ 2( £361 + £362 ) 

7t 

(b) For sin02 =± 1, 02 is a constant. For sin02 = 1, 02 = “■ However, if 

TC 

sin02 = -1, then 02 = For this case, if ( sin0isin02sin03 + cos03cos0i ) > 0, 
then 

01 = sin -1 2 ( £2 £3 + £i£ 4 ) j . (A. 6) 

If ( sin0isin02sin03 + cos03cos0i ) < 0,then 

01 = tc — sin -1 J 2 ( £ 2 E3 -+- £1 £4 ) j . (A.7) 

For this entire case, 03 = 0 . 
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Abstract 

We treat the question of large rotational maneuver and vibration 
stabilization of NASA Spacecraft Control Laboratory Experiment sys- 
tem (SCOLE). The mathematical model of SCOLE system includes 
the dynamical equations for rigid body slew maneuver and three- 
dimensional vibration of the rigid shuttle, the flexible beam and the 
reflector with an offset mass. The design approach taken here is to 
decompose the rigid mode control from vibration stabilization. Feed- 
back input (Shuttle torque)-output (attitude angles) map linearization 
technique is used for designing attitude control system for large-angle 
slewing. Linearization of input-output (i-o) map is accomplished by 
nonlinear inversion theory. It is shown that attitude control system 
assymptotically decouples the flexible dynamics and linear feedback 
law is easily designed for vibration suppression. For the synthesis of 
the control law an observer is designed. Simulation results are pre- 
sented to show in the closed-loop system large angle maneuver can be 
accomplished using only the measured state variables . 
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1. Introduction 


An important class of shuttle deployed payloads consists of cantilevered 
beam-like structures with massive tip bodies. A design challange was posed 
by Taylor and Balakrishnan [1] for the control of such a system. In the design 
challange, the Spacecraft Control Laboratory Experiment (SCOLE) was set 
up by NASA Langley Research Center to provide a standerd configuration 
to test control laws for investigators. 

The laboratory apparatus includes the Space Shuttle (orbiter) connected 
by a long (130ft) flexible beam to hexagonal antenna or reflector, as shown 
in Fig. 1. The reflector and Shuttle are treated as rigid bodies. When 
the Shuttle maneuvered, the flexible modes of the beam are excited. Two 
force actuators and a torquer have been provided to exert sufficient force to 
suppress the unwanted vibration of the beam. 

The design challange attracted the attention of many researchers. Re- 
search related to both dynamics and control system design has been per- 
formed and reported in litrature and SCOLE workshops at NASA Langley 
Research Center [2-10]. Nonlinear invertibility theory has been used in [11- 
12] for designing controller for a simplified spacecraft system. For such a 
system adaptive control system and sliding mode controllers have been also 
designed to compensate for parameter uncertainty and disturbance torques 
acting on the spacecraft [13-14]. 

We treat the question of large- angle rotational maneuver and elastic mode 
stabilization of SCOLE system based on input (Shuttle torque)-output (at- 
titude angles) feedback linearization for attitude control and linear feedback 
for vibration suppression. An i-o feedback linearizing control law is derived 
by using the inversion theory for nonlinear dynamical systems. The in- 
verse technique gives rise to decouple linear attitude dynamics and allows 
independent control of attitude angles. The use of inversion theory has an 
additional advantage since the elastic dynamics of the SCOLE configuration 
representing transverse vibration in two orthogonal plans and the torsional 
deformation of the elastic beam are assyptotically decoupled from the rigid 
mode dynamics. Since as usual only small elastic deformations are assumed 
here, the decoupled elastic dynamics are linear. Using linear control theory 
stabilizer is then designed for elastic mode stabilization. For the synthesis 
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of the controller only attitude angles, angular rates, tip elastic deflection 
components and torsion deflection at the tip of the beam are assumed to be 
measured by sensors. An observer is designed to estimate the elastic modes 
and their derivatives, and these are used for synthesizing the control law. 
Although the design approach has similarity to that of [11], this study differs 
in several ways. The mathematical model considerably differs from that of 
[11] since the end body combined here has mass and inertia also torsional 
modes are included. Furthermore, representation of elastic deformation dif- 
fers from [11]. The design of observer has led to realistic implimentation 
of control law in this paper using measured variables. The control of [11] 
requires knowledge of the complete state vector. 

The organization of the paper is as follows. Section 2 presents the math- 
ematical model. Inverse attitude control law is derived in section 3. Section 
4 presents elastic mode stabilizer. Observer design is presented in section 5 
and finally section 6 presents simulation results. 

2. Mathematical Model of SCOLE System 

We shall be interested in the rotational and elastic dynamics of SCOLE 
system in this study. To describe the orientation of the orbiter a body fixed 
orthogonal coordinate system with axes x,y,z is utilized. The attitude of the 
orbiter is defined by a sequence of rotation 0\,02, 0 3 (roll, pitch, yaw), with 
respect to an inertial coordinate system. The body-fixed coordinate frame 
has its origin at the point of attachment of the flexible appendage with the 
rigid Shuttle for this study. 

The orientation of the orbiter is completely described by the differential 
equations relating the rotation angles and the angular velocity components 
of the orbiter which are 

6 = M~ t u (1) 


where 


M~ t = (M r )- 1 


M t = 


cos02cos03 sin0 3 0 
—cos02sin0 3 cos0 3 0 
sin0 2 0 1 


and 
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The elastic beam is allowed to' undergo transverse elastic deformations along 
axes x and y, and torsional deformation about axis z. Elastic deformations 
are assumed to be linear combinations of admissible functions (f> x i(s), (f> y i(s) 
and i=l,2,3,...,n, where s denotes the distance on the beam from 

the attachment point of the beam and the orbiter. The chosen admissible 
functions are 


, / \ A • 71 Ct «' S ^ , Of,S 

<Pxi( s ) — AxiSin— + B x iCos-j- + CxiSinh— + D x iCosh— (2) 
L/ L/ Li JL 


Ct ’ S Qt'S Ot'S Ct’S 

<f>yi(s) = Ayisin—- -f Byicos-j- + C y isinh-~ + D y {Cosh (3) 
Li 1j Jj Li 


s s 

— A.ypxSXTiCt^i^r "f B^iCOSCC^i 

Li Li 


( 4 ) 


The parameters A x i, v4 y ,-, a, etc, are explained in [5]. 

The elastic deformations are explained as linear combinations of admis- 
sible functions as 

n 

( 5 ) 

( 6 ) 




»=1 

n 


«v( 5 ) = 

i=i 


u xp ( s ) — < ^'»'( 5 )?«(0 ( 7 ) 

1=1 

where u x , u y are the elastic deformations along x, y axes and u ^ is the tor- 
sional deformation about z-axis respectively. Here g,( t), i=l,2,3,...,n , are the 
modal deformation co-ordinates for the beam. Define u = (u>i,u; 2 ,u >3 ) T ,q = 
(?i>?2, ••M?n) r (T denotes transposition). Let Go = (Goi,Go 2 ,Gos) t be the 
external torquue vector applied on the orbiter, F = (F x ,F y ) T be the exter- 
nal control forces applied at the reflector center of the mass and M$ r be the 
moment applied about z-axis at the reflector center of mass. 

The equations of motion are obtained by using Lagrangian approach. For 
this, the expressions for the kinetic and potential energies are first obtained. 
The potential energy includes the strain energy due to transverse bending 
as well as due to torsional deformation of the elastic beam. The complete 
equations of motion describing rotational and elastic behavior of the SCOLE 
system have been derived in [5] (Readers may refer for the details). In this 
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study we shall be interested only in the rotational and elastic dynamics. 
These are 


7o A? 

tli 


+ 

' B x ‘ 

A% A 3 

. 9 ~ 

-Kq 

. . 


u 


(8) 


where 6 = ( 0 i, 02 5 03 ) r > K is the stiffness matrix, N 2 is a nonlinear vector 
function, B,- {i=l,2,3,...,n) are the constant input matrices, u = (Gj, F T , G 

R? is the control input vector and 

B \ = [/3 x 3 -#12] 

B2 = [O n x 3 -B22] 


Here I and O denote identity and null matrices of indicated dimensions. The 
matrix D 


D = 


Iq A2 
A\ A* 


is possitive definite symmetric matrix. 

Defining the state vector x=;(0 T ,ui T ,q T , q T ) T € R = (2 n + 6), one 
can obtain a state variable representation from (8) of the form 


x = f(x) + Bu 


where 


/(*) = ((Af ~ r a>) J , fZ(x), q J , fg (x)Y 


(9) 

( 10 ) 


f u (x), f q (x) and B are easily determined by comparing (8) and (9). 
Let the controlled output vector be 


y = 0 = c(x) 


( 11 ) 


and y c (t) € R 3 be any given reference trajectory. We are interested in deriv- 
ing a control law such that in the closed-loop system the output y(t) tracks 
the reference trajectory y c (<), and elastic modes are stabilized. This design 
objective will be accomplished by designing an inverse attitude controller for 
the orbiter for attitude angle tracking using orbiter moments; and an elastic 
mode stabilizer for vibration suppression using tip body external forces and 
moments. 
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3. Inverse Attitude Control 

In this section the design of control system for attitude trajectory tracking 
is considered. To derive the control law the inversion of the input (Go)- 
output ( 9 ) map is considered. In the following, an inversion algorithm [15- 
16] will be used to obtain an invertible system which will give the structure 
of the controller. This algorithm essentially gives a systematic procedure 
for obtaining a sequence of systems by the differentiation of the output and 
applying nonlinear transformations. Although, the algorithm is applicible 
for more general systems, the actuator dynamics have been neglected here 
for simplicity. 

Let B a be the first three columns of matrix B and B = [B g B r \. Then 
the system (9) can be written as 

X(t) = f(X,U r ) + BgUg (12) 


where 


u — (u 


T 
9 ’ 



Ug = Go 


= (F T , M^rY 


and 


/(X,lt r ) = f(x) + B r lt r 

In the following we define certain operators which will be useful in the 
sequel. For the vector function c(x), we define 

Lj(c)(x) = [dc(x)/dx]f(x) 

L}(c)(x) = L/(I,(c))(») (13) 

L Bg L f (c)(x) = [dL,(c)(x)/dx)B g 

Differentiating y and using (12) gives System 1 and System 2 of the form 
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System 1: 

x(t) = f(X,U T ) + BgUg 

Z i = Lj{c){x) (14) 

System 2: 

X(t) = f(l, U r ) + BgUg 

z 2 = L 2 j(c)(x, u r ) + D*(x)u g (15) 

where L 2 j(c)(x) is easily obtained using definition (13), z\ — y, z 2 — y, and 
D*(x) = LB g Lj(c)(x). Let R be the 3x3 submatrix formed by the first three 
rows and first three columns of D~ x . Then it can be easily seen that 

D* = M~ t R (16) 

We shall be interested in a region Cl of the state space where fi= {r. 
$2 ^ v ± tt/ 2}. We notice that R is nonsingular and M~ T exists in Cl. Thus 
IT is invertible in Cl and each component of the output has relative degree 
2. The inversion algorithm terminates here; and System 2 is invertible. 

Following [11], in view of (15), one obtains an i-o feedback linearizing 
control of the form 

Ug = d'-'[-l){c){x ) - p 2 e - p 1 e - p 0 x a + e c ] ( 17 ) 

where P, = diag(pij ), i=0,l,2; j=l,2,3\ 9=(9 — 9 C ) is the tracking error and 
x a is the integral of the tracking error, i.e, 

x, = 9 (18) 

The nonlinear control law includes proportional, derivative and integral 
(PID) type feedback. Integral feedback of trajectory error introduces robust- 
ness in the control system. 

Substituting control law (17) in the output equation (15), gives 


9 + P 2 9 + P\9 + Pqx s — 0 


(19) 
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which can be differentiated to yield 

6 +P 2 ~6 + pj + P 0 9 = 0 (20) 

It is obvious from (20) that each component 0,-(f)=(0,(<) — 0 tt -(t)) can be 
independently controlled where 0 C = (0 c i, 0c2, 8ds) T • Let the characteristic 
polynomial associated with 0, in (20) is of the form 

(s + A d )(s 2 + 2Cc.Wn«S +^L) = 0 ( 21 ) 

The matrices Pi are easily determined from (21). The parameters A^-, 
and ur nc i are chosen such that (20) is assymptotically stable. 

The inverse attitude control law causes tracking of the command trajec- 
tory 0 C . Indeed if matching conditions 

0 ( 0 ) = 0 

6 ( 0 ) = 0 (22) 
1 ( 0 ) = 0 

are satisfied; one has 0(t) = 0 c (t), t> 0, and 6 c (t) is exactly reproduced. 
However, slewing of the spacecraft excites the flexible modes and it becomes 
essential to damp the elastic oscillations so that the target in the space is 
tracked. 

4. Elastic Mode Stabilizer 


We assume in this section that the reference trajectory is such that 
0 c (t) — » 0*, where 6* corresponds to a desired terminal orientation of the 
orbiter and 0 C — »0, 0 C — >0, as t — > oo. For such a trajectory, in the closed-loop 
system, one has 6(i) — ► 6 C — ► 0*, and 0(<) — >0, u>(i) — >0, 0(t) — ► 0, as t — * oo. 

Interestingly, in the closed-loop system, in view of (8), the elastic dynamics 
are assymptotically decoupled, and reduces to 


A^q — —Kq 4- B^Ur 


Defining z=(q T , q T ), one obtains a state variable for (23) as 



Onxn 

L»xn 

z + 

0 

i = 

. -A^K 

Onxn 

y4 3 1 i?22 


(23) 

(24) 
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= A e z + B e u r 

For the stabilization of (24), one can use linear quadratic optimization or 
pole assignment technique. We have chosen pole assignment technique in 
this study. The closed-loop poles are chosen to be located at -.2±.lj, 
-.2±.3j, -.2±.5j, and -.2±.7j 
The linear feedback control law is 

u r = -F t z (25) 

The closed-loop system matrix A ec = ( A e — B e F e ) is a Hurwitz matrix. 

It is interesting to note that use of inverse attitude controller simplifies 
the stabilization problem to a linear design problem of a lower order system 
and this design is carried out independently. 


5. Observer Design 


The control laws (17) and (25) require knowledge of complete stste vari- 
ables. We assume that only measured variables are 6 , u, and the elas- 
tic deformation u x (L), u y (L ), and u^(L) at the tip of the beam. For the 
synthesis of the controller, it is essential to construct the variables q and 
q. With these measurements, the system (8) can be considered as an on- 
line linear time-varying system since the nonlinear function is of the form 
N 2 = N2c(0,9,u)) + N2j(9,6,u>)q which is linear in q and N 2 / depends on 
measured variables. However, design of observer for a time varying system 
is not easy. Since N 2 contains terms of second order in derivatives of 0,, and 
<7,-; N 2 / is small and , therefore, it can be neglected for simplicity in design 
of the observer. 

Setting N 2 jq= 0 in N 2 in (8); one obtains the following equation for elastic 
modes 

q = R2 

where R 2 consists of last n rows of D~ l . Since N 2c , and u are known, a state 
variable representation of (26) is 

Lk 0 2 + «(».»■“>.«) 




z = 


(27) 
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where g is obtained from (26), R 2 = [#21 R22], and R22 is n x n matrix. 
The associated measurement equation is linear function of q of the form 


Vm = c m q 


Define A m as 


Am — 


0 I 
—R22K 0 


then the state observer takes the form 


(28) 


i = A m z + g(S, 6, id, u) + H(y m - C m £) (29) 

where z is the estimate of state z, and the matrix H is determined such that 
the matrix (A m — HC m ) is Hurwitz. Again, the matrix H can be found using 
LQR or pole assignment technique. In this study, we use the pole assignment 
technique. 

The observation error z = (z — z) satisfies 

k = (A m -HC m )z (30) 

which is assymptotically stable and z(t) — » z, as t — » oo. For simulation 
matrix H is chosen such that the observer poles (eigen values of [A m — HC m ]) 
are at —8 ± 4 j, —8 ± 6;, —8 ± 8 j, —8 ± lOj. 

For the synthesis of the controller, q , q in control law (17) and (25) are 
replaced by q and q. The complete closed-loop system is shown in Fig. 2. 


6. Simulation Results 


We present here the results of digital simulation for various initial con- 
ditions and parameters. The appendix lists the maximum limits of control 
inputs system. The mass of the antenna is 400/6s, and the mass of the Shuttle 
body is 2.050 x10 s lbs [8]. 

For tracking a representative command reference trajetory is generated 
using a third order filter 

Oci +(2Cct^nct + A„)0 + 2(Cc»U>net A„ + ^nciWci + W nct' Act'^ci = ^nci'Act#* (31) 

such that its poles are at -Ac,- and {— Cct'Wnct ±jw nc i(l — P ara - 

maters chosen are = .707, and u >„„■ = A c/Cd- By a prorer choice of A CT -, 
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one can obtain desireable fast reference attitude trajestories for treking. The 
terminal angles 6* axe chosen 4°, 10°, and 25° (roll, pitch, and yaw respec- 
tively). 

(a) Trajectory Control: Stabilizer Open Loop 

In order to observe the behavior of the closed-loop system (8), 
(17), and (18), this simulation was done without the stabilizer. 

In this case the forces and moment M ^ are zero. The initial 
conditions were set as 0(O)=O, u>(0)=0, x 4 (0)=0, q(0)~ 0, ?(0)=(- • 
.l,-.l,-.l,-.l), q = q = 0. The elastic mode error (?(0) — ?(0)) was 
introduced to examine the state reconstruction ability of the ob- 
server. Only measured and estimated states are used for synthesis 
of the contrller. Selected response plots are shown in Fig. 3. The 
response time for the 6 is nearly 20 seconds. The simulation re- 
sults show the periodic and bounded oscillations of elastic modes 
and control moments as expected. 

(b) Trajectory Control: Stabilizer Loop Closed 

The complete closed-loop system including the stabilizer (24), and 
(25) was simulated. The initial conditions and inverse controller 
parameters of case (a) were retained. The selected simulation re- 
sults are shown in Fig. 4. The attitude angle response remained 
very similar to those of case (a). The stabilizer was switched at 
20 seconds. The stabilizer exerts a force in both x and y axes 
which causes a sudden change in responses when it is switched 
on. The poles of the stabilizer are properly selected such that 
the elastic oscillations are rapidly suppressed. Although, one can 
design a stabilizer for even faster stabilization, this will require 
larger control forces and moments. We notice only a small elas- 
tic deformation. Of course, the actual deformation will be more 
since the contribution of residual modes need to be added. We 
observed smooth attitude trajectory following and elastic mode 
stabilization. 


7. Conclusion 


The control of slewing maneuvers of NASA SCOLE system was consid- 
ered. The design was accomplished by decomposing the slew maneuver prob- 
lem from the elastic mode stabilization problem. Orbiter orientation control 
system was designed using nonlinear i-o map inversion theory. Using the 
inverse attitude controller any smooth attitude trajectory cam be followed. 
The inverse controller includes PID feedback of attitude errors. Stabiliza- 
tion of elastic modes was accomplished by a linear stabilizer using end body 
forces and moment. Interestingly, the design of stabilizer was carried out off- 
line seperately. For the synthesis of the controller, a linear state estimater 
was designed using elastic deformations of the beam at the tip. Extensive 
simulation results showed that the large maneuvers of the spacecraft can be 
performed to follow precise attitude trajectories and elastic modes can be 
stabilized using only output feedback. 
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9. appendix 


• CONTROL LIMITS 

Go{t) = 100,000 ft-lbs (moment applied at the orbitr mass center) 
Fo(t) = 800/6s (force applied at the orbiter mass center) 

M^r — 100,000 ft-lbs (moment applied at the reflector c.g .) 

• FIRST FOUR FLEXIBLE MODES OF SCOLE MODEL 


THREE DIMENSIONAL MODE SHAPE CHARACTERISTICS 


lu 


WHICH 


FREQ. (Hz.) I 0.27804240E+00 


0.120 12084E+01 
0.16282665E+00 
-0.19670286E+00 
-0.16983450E+00 
0.19616259E+00 
-0.10274618E-01 
0.57579 133E-02 
0.1 181005 7E-01 
-0.5 722046 2E-02 
0.19360955E-01 
-0.50748354E-01 
0.1 39780 18E-04 


MODE No. 




0.81300189E+00 


020540387E+01 
0.40868 188E-01 
-0.61958845E-01 
-0.41309992E-01 
0.61880796E-01 
-032438404E-01 
036509234E-01 
03439044 7E-01 
-0.364 64758E-01 
036611842E-O1 
0.9269890 IE-01 
-0.87320799E-05 


0.31357296E+00 


0.1275651 8E+01 
03885529 IE-02 
-0.14998387E-01 
-0.4332101 8E-02 
0.14985820E-01 
0.142 1978 1E+00 
-0326 9579 7E+00 
-0.19283 105E+00 
03264456 1E+00 
031835058E-01 
0.311 15282E-01 
-0.7599233 7E-05 


0.11856099E+01 


03480468 7E+01 
0.80641794E-01 
-0.67233377E-01 
-0.809 13938E-01 
0.671 063 16E-01 
0.13728679E+00 
-0.1 174693 2E+00 
-0.14085209E+00 
0.1 172505 7E+00 
0.82557693E-01 
-0.16158934E-03 
0.1043771 8E-07 
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• MATRICES A 3 AND A% 


A 3 = 


0.45879E+2 
0.36305E-1 
— 0.89042E— 1 
-0.14067E0 


0.363 05E-1 
0.621 IE +2 
0.11263E0 
-0.1471E0 


— 0.89042E— 1 
0.11263E0 
0.32737E+2 
— 0.6392E— 1 


-O.14067E0 
-0.1471E0 
— 0.6392E— 1 
0.2547E+3 


— 0.2133821E0 -0.3687057E+3 
0.380892 IE +3 -0.3030935E+2 
— 0.1808478E+3 -0.1318596E+3 
0.1423380E+3 -0. 113585 IE +1 


— G.7253901E— 1 
— 0.8427658E— 1 
-0.125799E0 
-0.236 7351E-1 


0.28202 17EO * 
* 24 “ 0.35746 92E0 
* 34 ® 0.2412807E1 
* 4 . 4 “ 0.52851 16E1 
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• CAPTIONS FOR FIGURES 
Fig. 1: 

SCOLE configuration showing the three-body model, coordinate sys- 
tem, and actuator control points. 

Fig. 2: 

Block diagram of the system. 

Fig. 3: 

Trajectory control: Stabilizer open loop 

(a) Attitude angles. 

(b) Control moments Gq. 

(c) Torsional deformations, u 

(d) State errors, (q t - fa), ( q 2 - fa). 

(e) Resultant tip defflection, ^Ju 2 + u 2 . 


Fig. 4: 

Trajectory control: Stabilizer closed loop 

(a) Control moments Go. 

(b) Control forces F. 

(c) Control moment M^ r . 

(d) Torsional deformations, u 

(e) Resultant tip defflection, ^u 2 -f u 2 . 
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and Vibration Damping of 
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and 
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Abstract 

The question of rotational maneuver and vibration stabilization of NASA 
Spacecraft Control Laboratory Experiment (SCOLE) System is considered. 
The mathematical model of SCOLE system includes the rigid body dynamics 
as well as the elastic dynamics representing transverse and torsional defor- 
mations of the elastic beam connecting the orbiter and end body (reflector). 
For the rotational maneuver, a new control law (orbiter control law) is de- 
rived using orbiter input torque vector. Using this control law detumbling 
and reorientation maneuvers of SCOLE system are accomplished, however, 
this excites the elastic modes of the beam. Interestingly the orbiter con- 
trol law asymptotically linearizes flexible dynamics. Using the linearized 
model, a linear feedback control law is designed for vibration suppression. 
An observer is designed for estimating the state variables using sensor out- 
puts which are used for the synthesis of the control law. Simulation results 
are presented to show that in the closed-loop system detumbling and re- 
orientation maneuvers can be accomplished and the effect of control and 
observation spillover is insignificant . 
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1. Introduction 

An important class of shuttle deployed payloads consists of cantilevered beam- 
like structures with massive tip bodies. A design challenge was posed by Taylor 
and Balakrishnan [1] for the control of such a system. In the design challenge, 
the Spacecraft Control Laboratory Experiment (SCOLE) was set up by NASA 
Langley Research Center to provide a standard configuration to test control laws 
for investigators. 

The laboratory apparatus includes the Space Shuttle (orbiter) connected by a 
long (130ft) flexible beam to a hexagonal antenna or reflector, as shown in Fig. 
1. The reflector and Shuttle are treated as rigid bodies. When the Shuttle is 
maneuvered, the flexible modes of the beam are excited. Two force actuators and 
a torquer have been provided to exert sufficient force and moment to suppress the 
unwanted vibration of the beam. 

Research related to both dynamics and control system design has been per- 
formed and reported in literature, NASA reports, and SCOLE workshops [2-10] 
at NASA Langley Research Center. Nonlinear invertibility theory has been used 
in [11-12] for designing controller for a simplified spacecraft system. Recently, this 
approach has been extended for the control of SCOLE system [12]. 

In this paper, Control System design for reorientation and detumbling maneu- 
vers and elastic mode stabilization of SCOLE system is considered. A new Orbiter 
Control Law is derived such that detumbling maneuver and attitude control of or- 
biter can be accomplished using shuttle input torque. Using this control law, any 
given reference angular velocity trajectory, w r , can be reproduced by the control 
system during detumbling maneuver. Furthermore, by a judicious choice of the 
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reference trajectory uv(t), desired orientation of orbiter can be attained inspite of 
the elastic oscillations of the beam caused due to maneuver of the space vehicle. 
In the closed- loop system including the orbiter control law, the elastic dynamics 
of the SCOLE configuration representing transverse vibration in two orthogonal 
planes and the torsional deformation of the elastic beam, are asymptotically de- 
coupled from the rigid mode dynamics. Using linear control theory, a stabilizer 
is designed for elastic mode stabilization. For the synthesis of the controller only 
attitude angles, angular rates, tip elastic deflection components and torsion deflec- 
tion at the tip of the beam are assumed to be measured by sensors. An observer 
is designed to estimate the elastic modes and their derivatives, and these are used 
for synthesizing the control law. 

The organization of the paper is as follows. Section 2 presents the mathematical 
model. The orbiter control law is derived in section 3. Section 4 and section 5 
present elastic mode stabilizer and observer design, respectively, and finally section 
6 presents simulation results. 

2. Mathematical Model of SCOLE System 

A body fixed orthogonal coordinate system with axes z, y, z is utilized to de- 
scribe the orientation of the orbiter. The body-fixed coordinate frame has its origin 
at the point of attachment of the flexible appendage with the rigid Shuttle for this 
study. The attitude of the orbiter is defined by a sequence of rotation 6 i,O 2 ,03 
(roll,pitch,yaw), with respect to an inertial coordinate system. 

The orientation is completely described by the differential equation 

9 = M~ t (6)u 

where T denotes transposition and 

« = («i i,u>2,u> 3 ) r , e = (<M 2 ,0 3 ) r , M~ t = (M t )~ x 


(1) 
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and 


M T , = 


cosQicosQz sin 0 3 0 

—cos$ 2 sin &3 cos $3 0 
sinQ 2 0 1 


Often the arguments of functions will be suppressed for simplicity. Equation (1) 


relates the rotation angles and the angular velocity components (u>i,a; 2 ,u> 3 ) of the 
orbiter. It is seen that M T becomes singular at 0 2 = ±(?r/ 2) for the chosen choice 


of rotations of the coordinate systems. However, this singularity can be avoided 


by a choice of different sequence of rotations of coordinate frames if it is required. 


In the following we shall be interested in the region of state space Ci in which 
$ 2 # ±*r/2. 


The elastic beam undergoes transverse elastic deformations along axes x and 
y , and torsional deformation about axis z. Elastic deformations are assumed to be 
linear combinations of admissible functions <£*;($), and ^,(s), i=l,2,3,...,n , 
where ‘s ’ denotes the distance on the beam from the attachment point of the 
beam and the orbiter. Here it is assumed that elastic deformations are adequately 
represented by n modes. The admissible functions <^ r ,(s), 4> y i(s) and <f>^i(s) are 


• / \ a • QjtS OtjS . - Qf|«§ _ i CK{$ /r»\ 

<t>xi{s) = A x {Sin— + tf X iCos— + C x i$inh— + D x iCOsh—- (2) 

L/ h L/ L 


i t \ a a ' S r> Q i S /-» • i a ' S rx i 0:,5 

0 y ,(s) = AyiSin-r- + Byicos— + CyiSinh-r- + Dyicosh—— 
Li Jj L/ Lf 

s s 

AipiSlTlQipi ” * 1 * COSQfj)i — 

Li Li 

The parameters A x {, A y i,ati etc, are explained in [6]. 


(3) 

(4) 


The elastic deformations are expressed as linear combinations of admissible 
functions as 

n 

“*(-») = <t>xi(s)qi(t) (5) 

t=l 
n 

u v( s ) =1 

1=1 


( 6 ) 
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«V-(s) = H <?W 5 )9«'(0 (7) 

i= 1 

where tt x ,u y are the elastic deformations along x,y axes and tty is the torsional 
deformation about z-axis, respectively. Here q,(t), i=l,2,3,...,n , are the modal 
deformation co-ordinates for the beam. Define q={q\,q 2 , ... ,q n ) T 1 the elastic modes 
vector, and Go = (C7 01 , £ 02 , Cos) 7 " the external torque vector applied on the orbiter. 
F = (F xy F y ) T is the external control force vector applied at the reflector center 
of the mass and M^ r is the moment applied about z-axis at the reflector center of 
mass. 


The complete equations of motion describing rotational and elastic behavior of 
the SCOLE system have been derived in [6] using Lagrangian approach (Readers 
may refer for the details). In this study we shall be interested only in the rotational 
and elastic dynamics. These are 


' lo 

a 2 

Cj 


Ur 

A 3 

. 9 . 



N 2 {0,6,u,q) 

-Kq 


+ 


B l 

B 2 


( 8 ) 


where K is the stiffness matrix, N 2 is a nonlinear vector function, B{ ( i=l,2 ) are 
the constant input matrices, u = (G % , F T , M^ r ) T € R? is the control input vector 
and 


B\ — [/3x3 # 12 ] 

B 2 = [O n x3 Bh] 


The matrices 7o, A, and B, have been derived in [6], and I and 0 denote identity 
and null matrices of indicated dimensions. The matrix D 


D 


Io A? 

A3 


is positive definite symmetric matrix. 



We are interested in deriving control laws such that in the closed-loop system 
detumbling and reorientation maneuvers can be performed. 

Detumbling Maneuver: 

Let a reference angular velocity trajectory u> r , be given such that u> r (t) — ► 0, as 
t — * oo. We are interested in deriving a control law such that in the closed-loop 
system u(t) = u>(t) — u> r (t ) — ► 0, as t — ► oo, and the elastic modes are stabilized. 
The derivation of orbiter control torque vector for u> r trajectory tracking will be 
based on nonlinear inversion of an input-output map. 

Reorientation Maneuver: 

Let 0 T {t) be a given reference attitude trajectory, such that 0 T (t) —* 0 m , a desired 
attitude angle, as t — » oo. We are interested in deriving a control law such that in 
the closed-loop system 9(t) —* 0 r (t), as t — ► oo, and elastic modes are stabilized. 
The rotational maneuver to attain desired orientation will be accomplished by a 
judicious choice of u T [t ) in the detumbling control law. 

3. Control Laws for Rotational Maneuver 

In this section, control laws for rotational maneuvers will be derived. 

3.1 Detumbling Maneuver: 

Here first the design of control system for u> r trajectory tracking is considered. 
To derive the control law the inversion of the input (Go)- output (u>) map is 
considered. For the inversion of input-output ( 1-0 ) map, we consider the derivative 
of u from (8) which is 

uj = I?ii[A r 2(0, 9, w, q) + B\u ] + R\-i[—Kq + B^u] (9) 

where D~ l =(Rij), i,j = 1,2 and Rij axe submatrices of D~ x of appropriate dimen- 
sions. Define u r =(F T , M^r) 7 , and u g = Gq. Then in view of the special form of 
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matrices i?,-, it follows that 

U> = 4) + B\2Ut} + R\2{ — Kq + #22Ur}] + R\\Ug (10) 

^ * 

= a m (6,0,u,q,q,u r ) + R u u g 

where 

a m = R\\{N 2 + Bnu r } + f?i2{~ K<1 + -# 22 * 4 } 

Since it is desired to track the reference trajectory u r , we choose a control law u g 
of the form 

u 9 — R\\ [~ a * _ Au>(t) -f uv(t)] (11) 

where A > 0. Substituting control law (11) in (10) gives 

t i>(<) = -\u(t) (12) 

Solving (12) gives, 

u}(t) = exp(— Xt)u>(0). 

Thus u(t) — ► u r (t), as 1 -too; and if w(0)=0, it follows that u>(t) = u> r (t), for all 
t > 0. By choosing suitable smooth trajectory u T (t) converging to zero, one can 
accomplish desirable detumbling maneuver. If u> r (t) — ► 0 as t — * oo; then 0(t) — > 0 
and 8(t) — > 0j, a constant vector. In such a case when detumbling maneuver is 
completed, the orbiter attains a fixed orientation. 

3.2 Reorientation Maneuver: 

Now we shall derive a control law u g such that 6(t) — ► 0 T (t), a given reference 
attitude trajectory, as t — + oo, giving a desired orientation of the orbiter. This will 
be accomplished by a suitable choice of reference trajectory uy(<) in the detumbling 
control law (11), such that the desired orientation is asymptotically attained. 
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Since in the closed-loop system (8) and (11), u;(t) — ♦ u> r (t), we have that 
— ►M -1 (0)u; r (t), as t — > oo. Thus the asymptotic relationship between 
9 and u> obtained from (1) is 

9 = M~ T (9)u r (13) 

as t — * oo. Since it is desired to track 0 r (t), in view of (13), it follows that a 
suitable choice of u r is 

uv = M T (e)[-tid + e T \ (i4) 

where [x > 0, and 9 = (6 — 9 r ). 

From equation (1), we have 

0 = M~ T {9){u + u r ) 

= M~ t (9)[u + M T (9)(-fi9 + 0 r )] (15) 

which yields 

'0 = + M- t {0)u (16) 

Since for any choice of u> T (t ), in the closed-loop system including the control law 
(11), w — » 0 as t — * oo, asymptotically (16) reduces to 

9 = -n9 (17) 

It is interesting to note that by a proper choice of reference trajectory u T as given 
in (14), one obtains a first order linear differential equation asymptotically for the 
attitude angle tracking error. It is obvious from (17) that for any w(0), w r (0), 0(0), 
and 9 r ( 0); 9(t) — ► 0, as t — ♦ oo. Let 9 — 9 * corresponds to a desired orientation 
of the orbiter. Thus if 0 T {t) — * 0*, then 9(t) —* 9*, as t — * oo, and the desired 
orientation of the orbiter is attained. 

Furthermore, with the choice of <D(0) = 0 and 0(0) = 0, it follows from (12) 

that w(<) = 0 for all t > 0. Then substituting u)(f )=0 in (16), and using 0(O)=O, it 
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follows that 0(t) = 0 for all t > 0, and in the closed-loop system trajectories w T (t) 
and 0 T (t) are exactly reproduced for all t > 0. This way, the SCOLE system can 
be maneuvered to attain specified time- varying orientation at each instant t > 0, 
which is useful for tracking a rotating target in space. 

The control law for reorientation maneuver is obtained by combining (11) and 
(14) which is 

“» = Rnl-«'-Hu-M T (e)(-rf + 6 r )} + M T {6,6)(-n6+e,) + M T (0){-n'd + 0,)] 

(18) 

where M(0,0) = [(dM,j/d0)0], i,j=l,2,3 and Mi, are the elements of the matrix 
M. 

We notice from (18) that the second order derivative of 0 T is required in the 
control law. For this purpose one can choose a third order attitude angle command 
generator of the form (i= 1,2,3) 

0ri +(2CriU> nrt - + Ari)^ri + 2(( r i^>nri Ki + ^nriW^i + W nn^r«^ r . = ^nri^rt#* (19) 

where 0 T = (#ri, 0r2> QriY • The parameters A r ,-, and u ; nri - are appropriately 
chosen to obtain desired reference trajectory 0 T {t) terminating at 6’. It can be 
easily varified that the poles associated with (19) are at -A r , and {— Cr«w n r« ± 

J^nr,(l - Cr») 1/2 }* 

4. Elastic Mode Stabilizer 

We note that when the detumbling controller (11) is used, u>(t) — + 0, and 
w(<) -* 0, as t — ► oo for an approriate choice of u > r (t). Also in the closed-loop 
system including the reorientation controller (18), one can choose 0 r (t ) such that 
0(t) — * 0 *, 0(t) —> 0, u(t) —> 0, and w(f) — > 0, as t — » oo. However, rotational 
majieuvers cause transverse and torsional vibrations of the beam. Thus it becomes 
necessary to design a stabilizer so that the elastic modes can be stabilized. 
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In this section the design of an elastic mode stabilizer is considered. Since 
in the closed-loop system (8) and (11) or (8) and (18), w(t) —* 0 and w(t) — ► 0, 
irrespective of the elastic oscillation of the beam; it follows from (8) that the elastic 
dynamics get asymptotically decoupled from the rigid dynamics as t — ► oo. Setting 
u>{t) = 0 in (8), gives the decoupled elastic dynamics of the form 


A^q = —Kq + B 22 U r 


( 20 ) 


For the damping of the vibration, it is necessary to stabilize system (20). 


Define z=(q T ,q T ), to obtain a state variable form from (20) as 


Onxn 

Inxn 

2 + 

0 

f 

1 

Ju 

W | 
£ 

Onxn 

^3^22 


( 21 ) 


= A t z + B t u r 

where matrices A e , and B e are defined by (21). For the stabilization of (21), linear 
quadratic optimization or pole assignment technique can be used. We have chosen 
pole assignment technique in this study. The linear feedback control law is 


u r = — F e z (22) 

The matrix F e is chosen such that the closed-loop matrix A ec = ( A e — B e F e ) is a 
Hurwitz matrix and has specified eigenvalues. It is interesting to note that the use 
of orbiter attitude controller simplifies the stabilization problem to a linear design 
problem of a lower order system and this design is carried out independently. 

Although, the control laws u g and u r have been derived independently and the 
derivation of u T is based on a reduced order linear model, boundedness of elastic 
modes can be established in the complete full order closed-loop system for any 
bounded trajectory u r (t) and u> r (t). In the following we shall obtain a bound on 
the elastic modes. 
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In the closed-loop system with control law u g of (11), we note that u — + 0 as 
i — ► oo. The elastic dynamics obtained from the systems (8), (11), and (22) are 
given by 

i = A ec z 4- B a u (23) 

where f?j=[0 (-Aj 1 ^ ) T ] T . Using u = u r — Xu obtained from (12), (23) gives 

i = A ec z + B a (u r - Xu) (24) 

The solution of (24) is 

z(t) = i i>(t - t 0 )z 0 + f <f>(t- r)B a (u r - X u)dr (25) 

J to 

where z 0 = z{to) and the transition matrix <f>(t) associated with the matrix A ec is 
given by 

4>(t) = exp(A ec t) 

We assume in the following that u T and u T are bounded. Furthermore, the 
solution of (12) gives, ||tD|| = e -A< ||u;(0)|| < oo, where ||.|| denotes the Euclidean 
norm of a vector. Thus it follows that there exists a function 7 (t) < 00 such that 
for all t > 0. 

||S 3 (u; r - Aw) || < 7 (t) < 70 < 00 

where 70 is a constant. Since A ec is a Hurwitz matrix \\<j>(t — <o)||i < me' 0 * 1 '* 0 ) 
for some m > 0 and a > 0, where ||.||j is the induced matrix norm [13]. Let 
Ud = u r — Xu and z(to) = zq . Taking the norm and using triangle inequality, (25) 
gives (*o > 0) 

B*(0ll < ll^(« - <o)IMI*bll + [' \m - r)|| i .||^(r)||dr 

Jto 

< me- o(t - to )||z 0 || Am f e~ a(t ~ r h (r)dr 

Jto 
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< me- 0 ^-‘°)|| 2 o|| + (m/a)( 1 - e^^bo 

< m(j| 2 0 || + Q _1 7 o) < oo (26) 

This proves that in the closed-loop system q(t) and q(t ) are bounded provided 
that the reference trajectory is such that u T {t) and i l> r (t) are bounded functions 
and (26) gives a bound on ||z(t)||. 

For a choice of a reference trajectory u r (t) in the closed-loop system with the 
detumbling control law, for which u> r (t) — > 0 and u r (t) — ► 0, as t — * oo, iu(t) and 
cj(t) — » 0 as well, as t — > oo. In such a case, in the closed-loop system, it is obvious 
from (23) that z(t) — ► 0 as t — ♦ oo, since A ec is Hurwitz matrix. 

The stabilizer has been designed by setting d>(t) = 0 in (8). Since u(t) and tu(f) 
asymptotically tend to zero for a proper choice of u> r (t), a question arises: at what 
instant should the elastic mode stabilizer-loop be closed? According to (26), if the 
stabilizer is switched at the instant to, following a rotational maneuver command 
beginning at t = 0, the elastic modes remain bounded for all t > 1 0 , provided 
z(t 0 ) is bounded. However, since the stabilizer is not in the loop during the period 
[0,< o ], divergent elastic oscillations will result during this period in the system 
(23) in which A ec =A e (since F e = 0, z’.e; zero stabilization signal) if the reference 
trajectories v T (t) happen to contain sinusoidal signals of frequencies coinciding 
with the natural frequencies of vibration of the beam. These natural frequencies 
of oscillation are given by the purely imaginary eigenvalues of the open-loop matrix 
A, e- Thus it is preferable to keep the stabilizer loop closed from the instant t = 0 
to avoid structural resonance. 

5. Observer Design 

In this section, the design of an observer to estimate the states is considered. 
The estimation of state is necessary so that the controller can be synthesized using 
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measured signals. In this study it is assumed that the only measured variables 
are 9 , u>, and the elastic deformations u x (L ), u y (L), and u^,(L) at the tip of the 
beam. The variable 9 is computed using (1). For the synthesis of the controller, 
it is essential to estimate the variables q and q. With these measurements, the 
system (8) can be considered as an on-line linear time-varying system since the 
nonlinear function is of the form 7V 2 = -/V 2c (0, 9, u) + iV 2 /(0, 9, w)g which is linear in 
q and Nij depends on measured variables. It is interesting to note that N? contains 
terms of second order in derivatives of 0,, and <?,; therefore, iV 2 / can be neglected 
for simplicity in design of the observer. We shall make this simplification, since 
designing observer for the nonlinear model is extremely difficult. 


Setting iV 2 /<7= 0 in iV 2 ; one obtains from (8) the following equation for the 
elastic modes 

^2 c(^i^w) + B\U 
—Kq + B 2 U 

where f? 2 consists of last n rows of D~ 1 . Since iV 2c , and u are known, a state 
variable representation of (27) is 


q = Ri 


(27) 


z = 


0 I 

— i? 22 /f 0 


z + g(9,9,u,u) 


(28) 


where g is easily obtained from (27), i? 2 = [i? 2 i /? 2 2 ], and /? 22 is an n x n matrix. 

The associated measurement equation is linear function of q of the form 


Vm ~~ Cm? 


(29) 


where y m =(u x (L), Uy{L), u^,(L)) T . C m is computed using (2)-(7). Since the system 
(28) is linear and the function g can be measured on-line, one can design an 
observer using linear control theory. 

Define A m as 


A 


m 


0 7 

— 7? 22 7f 0 
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then the state observer takes the form 

z = A„iZ + g(O,0,u;,v,) + H(y m - C m z) (30) 

where z — ( q T , q ) T is the estimate of state z , and the matrix H is determined such 
that the matrix (A m — HC m ) is Hurwitz. In this study, we use the pole assignment 
technique to obtain feedback matrix H. 

The observation error z = (z — z) satisfies 

'z = (A m - HC m )z (31) 

which is asymptotically stable and z(t) — + 0, as t — ► oo. 

For the synthesis of the controller, q, q in control law (11), (18), and (22) are 
replaced by the estimated variables q and q. The complete closed-loop system is 
shown in Fig. 2. 

6. Simulation Results 

We present here the results of digital simulation. The numerical values of 
various matrices /o, A 1, K, J5j, B 2 and C m , and the complete expression of 
nonlinear function -/V 2 have been taken from [6]. Here we set n = 8. The appendix 
lists the maximum limits of control inputs. The mass of the antenna is 400/6s, and 
the mass of the Shuttle body is 2.050 x10 s /6s [8j. 

We chose A = ft = 0.55 for reorientation control law and detumbling control law. 
For tracking, a representative command reference attitude trajectory is generated 
using filter (19) with parameters A r , = 0.55 (r. = .707, and u? nr , = A „■/£„•. The 
terminal angles 0 * ( i=l,2,3 ), are chosen as 15°, 14°, and 13° (roll, pitch, and yaw, 
respectively). Since in the practical situation, the control magnitudes are limited, 
we have introduced control saturation functions in digital simulation. In this study, 
the initial conditions are 0(0) = 0, and g(0) = g(0) = 0. The observer poles are 
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placed at -18 ±;'1.2424, -18 ±>2.1571, -18±>4.7615, -18 ±>7.4541, -18 ± 
>12.5634, -18 ± >29.9435, -18 ± >34.3916, and -18 ± >77.1321. The stabilizer 
poles are located at —1.2 ± j’ 1 .6841 , —1.2 ± j 1 .8951 , —1.2 ± >4.8035, —1.2 ± 
>7.4018, -1.2±>12.4961, -1.2±>29.9426, -1.2±>34.3913, and -1.2±>77.1313, 
by using pole placement technique. These pole locations have been selected by 
trial and error and by observing the simulated response. 

(a) A Detumbling Maneuver: 

For detumbling maneuver, an exponentially decaying reference trajec- 
tory of the form u> T (t) = u> T (0)exp(— 0.5t) for tracking was consid- 
ered. The initial conditions were u>(0) = (3, 2, 1) T degrees/sec, and 
w r (0) = (5,4.5,4) r degrees/sec. The initial conditions of the observer 
were 9,(0) = (—0.01) and 9,(0) = 0, i=l,2,....,8. Thus g(0) = (0.01) 
and initial tracking error w(0) = (— 2, — 4.5, — 3) T degrees/sec. Selected 
responses are shown in Fig. 3. It is interesting to observe that inspite 
of the control saturation and the use of only estimated states in the 
control law, the tracking error tends to zero in about 10 seconds. We 
observe that 9 converges to a constant value as predicted, and elastic 
modes oscillations are suppressed. 

(b) A Reorientation Maneuver: 

Simulation was done to examine the reorientation maneuver capability 
of the controller. For this purpose, the tracking of the reference attitude 
trajectory, 6 r {t), terminating at 8 m (t) = (15°, 14°, 13°) T was considered. 
The initial conditions were the same as described in case (a). The 
selected responses are shown in Fig. 4. The response of an estimated 
state was somewhat similar to those of case (a). Here the attitude angle 
response 9(t) closely followed the reference attitude trajectory 0 T {t), but 
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state was somewhat similar to those of case (a). Here the attitude angle 
response Q(t) closely followed the reference attitude trajectory 0 r (t), but 
due to the control limits and the state estimation error, there was a small 
error in 6 tracking. The attitude angle 6 converged to 6" and elastic 
oscillation was suppressed in less than 10 seconds. 

(c) Effect of the Modal Error: The Spillover: 

Simulation was done to examine the effect of control and observation 
spillover by assuming that the actual dynamical model of SCOLE Sys- 
tem has 10 elastic modes qi(t). However, the controller, designed for the 
lower order model having only 8 elastic modes, was retained. To include 
the effect of observation spillover, the measured signal y m (t ) in the ob- 
server (30) was replaced by y m = C ma q a , where q a = (q u q 2 , ...., qio) T ■ 
Here C ma is a 3 x 10 output matrix whereas C m is a 3 x 8 matrix. 
The initial conditions of case ( b ) were retained except 9,(0) = (—1.1), 
i=l,2,3,....,8; u>(0) = 0; and 9,(0) = 9 = 0, i=l,2,3,....,10. Selected 
responses are shown in Fig. 5. We observe only a small effect of model 
truncation on on the controller. Stable responses were obtained. The 
desired orientation is attained in about 30 seconds. Bounded and con- 
vergent oscillation of the beam was observed in the simulation. These 
oscillatory responses are insignificant after 100 seconds. However it 
should be noted that the structural damping (which has been neglected 
in the study) will cause the dissipation of energy and eventually even 
these insignificantly small oscillations will vanish. 

7. Conclusion 


The control of slewing maneuvers of NASA SCOLE system was considered. 
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The design was accomplished by decomposing the slew maneuver problem from 
the elastic mode stabilization problem. Orbiter control system were for detum- 
bling and reorientation maneuvers. The orbiter attitude controller was obtained 
by a judicious choice of reference angular velocity trajectory in the detumbling 
control law such that any smooth attitude trajectory can be followed. Stabiliza- 
tion of elastic modes was accomplished by a linear stabilizer using end body forces 
and moment. For the synthesis of the controller, a linear state estimater was de- 
signed using elastic deformations of the beam at the tip. Extensive simulation 
results showed that rotational maneuvers of the spacecraft can be performed to 
follow precise angular velocity or attitude trajectories and elastic modes can be 
stabilized using only output feedback. 
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9. APPENDIX 


• CONTROL LIMITS 

G 0 {i) = 100,000 ft-lbs (moment applied at the orbiter mass center) 
Fo{t) = 800/6s (force applied at the orbiter mass center) 

M^r — 100,000 ft-lbs (moment applied at the reflector c.g.) 

• CAPTIONS FOR FIGURES 


Fig. 1: 

SCOLE configuration showing the three-body model, coordinate system, and 
actuator control points. 

Fig. 2: 

Block diagram of the system. 

Fig. 3: 

Detumbling Control Law: 

(a) Attitude angles. 

(b) Control moments Go. 

(c) Control forces F. 

(d) Angular velocity u>. 

(e) Tracking error u> 

(f) Control moment M^ r . 
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Fig- 4: 

Reorientation Control Law 

^ Attitude angles. 
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Fig. 3(a) 
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Fig. 3(b) 
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Fig. 3(c) 
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Fig. 3(d) 







- 117 - 



Fig. 3(e) 
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Fig. 3(f) 
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Fig. 3(g) 
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Fig. 3(i) 
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Fig. 4(b) 
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Fig. 4(c) 
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Fig. 4(d) 





Moment, Reflector C.G., ft-lbsxlO 3 


126 



Time, Seconds 


Fig. 4(e) 
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Fig. 4(g) 
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Fig. 5(a) 
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Fig. 5(b) 
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